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FORWARD 

The  research  reported  herein  was  performed  by  HYDRONAUTICS, 
Incorporated  for  Cornell  Aeronautical  Laboratory,  Inc.  (CAL), 
and  supports  the  Soil  Mechanics  Task  of  the  Off-Road  Mobility 
Research  (ORMR)  program  sponsored  by  OSD/ARPA  for  which  CAL  is 
the  prime  contractor  and  the  U.  S.  Army  Research  Office  -  Durham 
is  the  contracting  agency. 

The  reported  study  is  part  of  a  broad  research  program 
designed  to  acquire  off-road  mobility  knowledge;  develop  analysis, 
prediction  and  decision  methodologies;  and  organize  knowledge 
and  methods  to  facilitate  use  by  military  planners,  vehicle  de¬ 
signers  and  field  personnel  concerned  with  off- road  mobility. 

Under  these  broad  objectives,  soil  mechanics  research  is  directed 
at  developing  techniques  for  analyzing  off-road  mobility  problems 
related  to  soft-soil  conditions.  The  present  study  is  concerned 
with  steady  state  two-dimensional  flow  generated  by  a  driven 
rigid  cylindrical  wheel  moving  in  soft  saturated  clay,  assumed 
to  behave  like  a  rigid  plastic  material.  The  study  was  per¬ 
formed  during  the  period  April  1968  -  September  1968.  Funds 
for  this  study  were  allocated  to  HYDRONAUTICS,  Incorporated  by 
CAL  under  ARPA  Order  No.  841  dated  7  May  1966  with  amendment 
Nos.  1-4.  Dr.  Patrick  Miller  was  the  CAL  ORMR  Soil  Mechanics 
Task  Leader;  Mr.  G.  E.  Bartlett  was  ORMR  Program  Manager  at  CAL; 
and  Mr.  A.  N.  Tedesco  was  ARPA  Technical  monitor. 
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Dotted  symbols  have  dimensions.  Undotted  variables 
are  dimensionless:  length  by  dividing  with  r  ' ,  velocities  by  ' , 
stresses  by  k'  and  forces  by  k'r*. 

g  Gravity  acceleration 

G  =  Dimensionless  gravity  plastic  number 

K 

h'  Thickness  of  the  plastic  zone  at  the  bottom  of 

the  wheel  (h  =  h'/tr') 

P'V* 

H  =  — — —  Dimensionless  Inertial  plastic  number;  also  hori¬ 

zontal  force  acting  on  the  wheel  (H  =  H'/r'k') 

J  Bessel  function  of  zero  order 

o 

k'  Yield  stress  (max.  shear  stress) 

t,L  Auxiliary  lengths  defined  in  Figures  17  and  19 

M'  Torque  on  the  wheel  (M  =  M  /k'r,a) 

n  Direction  normal  to  streamline:; 

p!  Isotropic  pressure  (p  =  p'/k1) 

r'  Wheel  radius 

r  Radial  coordinate 

R'  Resistance  (R  =  F'/k'r  ) 

s  Direction  along  streamlines 

V  !  -  U  '  \ 

S  Slippage  coefficient  S  =  —  -  ,  | 

,  w  / 

u',v'  Velocity  components  in  cartesian  coordinates 

(u  =  u7Uo',  v  =  v’/U0'  ) 

u 

b 


Velocity  in  the  plastic  zore  at  the  bottom  of 
the  wr.eel 
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li  '  Translational  velocity  of  the  wheel 

o 


V' 


V  '  =  uu'r' 
w 


Va’V/3 

Vve 

W' 


x'  ,y' 


X  |X„ 

a  0 

z ' 

a,0 

7' 

6V,6T 

xx  yy  'xy 
• 

e 

x 

p' 
a 1 

a  ,  a 

x  y 


w 


Velocity  modulus  (V  =  V'/U  1 ) 

Peripheral  rotational  velocity  of  the  wheel 

(V  =  V  '  /U  '  ) 
v  w  w  o  ' 

Velocity  components  in  characteristic  directions 

Velocity  components  in  radial  coordinates 

Vertical  force  acting  on  the  wheel  (W  =  W'/lc'r') 
Cartesian  coordinates  (x  =  x'/r'»  y  =  y'/r') 
Curvilinear  characteristic  coordinates 

Slnkage  (z  =  z'/r') 

Characteristic  directions 
Shearing  strain 

Velocity  and  shear  stress  jumps 
Rates  of  strain 

Auxiliary  angle  defined  in  Figure  19 

Angle  between  streamline  and  x  direction  (in  par- 
tlcu1ar  polar  coordinate  along  the  wheel) 

Slope  of  characteristic  directions 

Soil  density 

Normal  stress  (a  =  a' A') 

Cartesian  components  of  the  normal  stress  tensor 
Normal  stress  acting  on  the  wheel 


CijOa  Principal  normal  stresses 

T'  Shear  stress  (t  =  T'/k') 
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T  Constant  shear  stress  along  the  wheel 

w 

0  Angle  between  a  characteristic  and  x  axis 

0  0  along  the  wheel 

w 

u) '  Rotational  speed  of  the  wheel 
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ABS  TRACT 

The  steady  two-dimensional  flow  beneath  a  driven  rigid 
cylindrical  wheel  moving  in  a  soft  saturated  clay  is  solved  by 
assuming  that  the  soil  behaves  like  a  rigid-plastic  material. 

The  general  equations  of  plastic  flow  are  discussed,  with  empha¬ 
sis  on  the  inertial  effects  which  generally  are  not  negligible. 

A  plastic  flow  pattern  is  suggested  and  the  quasi-static  equa¬ 
tions  of  flow  are  integrated  in  two  regions  of  the  plastic  zone. 
An  approximate  solution  provides  the  magnitude  of  the  recovery 
angle,  the  vertical  and  horizontal  forces  and  the  torque  acting 
on  the  wheel  as  function  of  wheel  radius,  sinkage  and  shear 
stress  along  the  wheel  (assumed  constant  on  the  bow  portion). 

A  minimum  slippage  necessary  to  maintain  the  shear  stress  on 
the  wheel  is  found.  The  theoretical  results  are  compared  with 
some  existing  measurements  and  the  agreement  is  generally 
satisfactory . 
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1.  INTRODUCTION 

The  capability  of  a  vehicle  to  navigate  in  a  soft  soil  de¬ 
pends  on  the  magnitude  of  the  forces  developed  at  the  wheel 
(or  track)  and  soil  interface.  The  prediction  of  these  forces 
as  functions  of  the  wheel  geometry,  soil  characteristics  and 
wheel  speed  is  necessary  in  order  to  design  vehicles  or  to  esti¬ 
mate  the  mobility  of  a  given  vehicle  under  different  conditions. 

Unfortunately,  at  this  stage  of  development  of  the  art  (see 
Section  2)  there  is  no  reliable  theory  which  can  produce  design 
formulae.  Moreover,  even  the  fundamentals  are  not  well  under¬ 
stood  and  similitude  criteria  which  allow  extrapolation  of 
model  tests  to  prototype  are  not  clearly  formulated. 

The  purpose  of  this  study  is  to  solve  In  a  basic  way  the 
problem  of  soil  flow  beneath  a  wheel  under  some  simplifying 
assumptions.  While  existing  approaches  use  empirical  results 
obtained  from  static  tests  of  plates  on  the  soil  surface  and  ap¬ 
ply  them  to  the  moving  wheel,  in  the  present  study  the  whole 
region  of  soil  flow  beneath  the  wheel  is  considered.  The  solu¬ 
tion  of  the  problem  follows  the  lines  of  the  classical  applied 
mechanics:  the  soil  behavior  is  represented  by  material  consti¬ 
tutive  equations,  the  stresses  and  velocity  fields  are  inter¬ 
related  through  the  equations  of  motion  and  the  particular 
solution  for  the  wheel  is  obtained  by  integration  with  the  ap¬ 
propriate  boundary  conditions. 
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Unfortunately,  the  difficulties  encountered  in  such  a  type 
of  solution  are  formidable.  The  most  important  of  them  is  the 
lack  of  a  good  mathematical  representation  of  the  soil  behavior 
in  form  of  constitutive  equations,  i.e.,  relationships  between 
stresses  and  deformations.  But  even  if  some  simplified  consti¬ 
tutive  equations  are  adopted,  the  exact  Integration  of  the  equa¬ 
tions  of  motion  in  the  wheel  case  seems  impossible  at  this  stage. 

For  the  above  reasons,  in  this  preliminary  study  the  problem 
has  been  simplified  by  adopting  the  following  assumptions:  (i) 
the  soil  behaves  like  a  rigid-plastic  incompressible  material, 
(ii)  the  wheel  is  two-dimensional  and  rigid,  (ill)  the  motion 
is  steady  and  (iv)  the  soil  surface  is  plane  and  unperturbed 
far  from  the  wheel.  Assumption  (i)  limits  the  applicability  of 
the  results  to  soft  clays  with  a  high  water  content,  although 
even  in  this  case  the  plastic-rigid  model  is  a  foregoing  simpli¬ 
fication.  It  may  describe  roughly  the  behavior  of  a  saturated 
undrained  clay  with  a  dispersed  structure  (Reference  l4,  p.  369) > 
(Reference  7,  p.  188). 

Considering  the  highly  empirical  character  of  existing 
theories  it  is  felt  that  the  present  work,  in  spite  of  the  above 
simplifications,  constitutes  an  important  step  forward  in  the 
understanding  of  soil-vehicle  interaction.  This  is  only  a  first 
step,  and  the  results  obtained  so  far  encourage  the  contin^ 
uatlon  of  the  study  along  the  same  lines,  with  steady  improve¬ 
ment  and  extension  of  results  to  more  complex  conditions,  like: 
representation  of  soil  as  a  work-hardening  material  and  incorpo¬ 
ration  of  friction,  consideration  of  three-dimensional  wheels  and 
of  unsteady  motions.  Influence  of  wheel  or  track  shape,  etc. 


HYDRONAUTICS,  Incorporated 


-4- 

2.  BRIEF  LITERATURE  REVIEW 

The  purpose  of  this  section  is  not  to  give  a  comprehensive 
review  of  the  publications  in  the  field  of  vehicle-soil  inter¬ 
action,  but  a  critical  discussion  of  the  most  important  contri¬ 
butions,  relevant  to  this  study. 

A .  Theoretical  Works  on  Soil  Wheel-Interaction 

The  most  known  and  applied  theories  of  wheel-soil  interac¬ 
tion  are  those  of  Bekker  and  Uffelraann. 

Bekker's  theory  (2),  (3)  is  based  on  the  following  assump¬ 
tions  on  development  of  stresses  on  the  soil-wheel  interface: 

(l)  The  normal  stress  at  a  point  is  related  to  the  wheel  sinkage 
at  the  same  point.  The  relationship,  known  as  Bernstein  equa¬ 
tion,  is  obtained  empirically  by  pressing  plates  in  the  soil  and 
representing  the  dependence  of  force  on  sinkage.  (li)  The  hori¬ 
zontal  resisting  force  to  the  wheel  motion  is  a  result  of  the 
work  done  in  compacting  the  soil  and  creating  the  rut.  There 
is  no  soil  recovery  beyond  the  wheel  bottom.  (ill)  The  shear 
stress  along  the  wheel  is  computed  by  using  the  relationship  be¬ 
tween  shear  stress  and  deformation  obtained  empirically  by  pulling 
a  flat  plate  on  the  soil  surface. 

The  soil  layer  adjacent  to  the  wheel  has  an  angular  defor¬ 
mation  which  is  computed  by  assuming  that  on  the  wheel  side  the 
speed  is  equal  to  the  wheel  rotational  speed  and  or  the  soil 
side  equal  to  the  translational  speed.  This  permits  the  deter¬ 
mination  of  the  shear  stress  as  function  of  slippage  and  position. 
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By  using  seven  soil  constants  in  the  representation  of  the 
empirical  curves  and  some  additional  computational  simplifica¬ 
tions  Bekker  was  able  to  calculate  the  forces  acting  on  the 
wheel  by  using  the  above  assumptions. 

The  inconsistency  of  Bekker' s  theory  when  compared  with  ex¬ 
perimental  results  have  been  largely  discussed  in  literature 
(^ee,  for  instance,  (19),  (20))  and  will  not  be  repeated  here. 

From  a  basic  point  of  view,  Bekker 's  theory  is  an  attempt 
to  apply  results  of  measurements  of  forces  on  flat  plates  in 
static  conditions  to  moving  wheels.  There  is  no  a  priori 
reason  that  such  an  attempt  should  be  successful  since  the  soil- 
flow  beneath  a  flat-plate  is  different  from  that  beneath  a 
wheel.  If  we  adopt  the  plastic -rigid  model  and  observe  the 
plastic  flow  pattern  beneath  a  footing  (7),  (8)  and  that  be¬ 
neath  a  rolling  surface  (l),  (8)  we  will  find  that  they  are 
different.  Moreover,  the  normal  stress  at  a  certain  point  is 
depending  on  the  whole  flow-field  and  not  on  the  vertical  dis¬ 
placement  of  that  particular  point.  There  is  no  theoretical 
ground  for  the  assumption  of  pressure-slnkage  relationship  along 
the  wheel. 

The  rolling  resistance  of  a  wheel  cannot  be  attributed  to 
the  soil  compaction,  excepting  possible  agricultural  soils 
with  high  content  of  air-filled  pores.  A  soft  clay  with  a  high 
water  content  is  practically  Incompressible  and  in  a  two-dimen¬ 
sional  case  there  is  no  ultimate  rut.  There  is  definite  evidence 
of  soil  recovery  at  the  rear  with  normal  stresses  contributing 
to  draw-bar  pull  (18),  (22). 
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The  distribution  of  shear  stresses  along  the  wheel  is  dif¬ 
ferent  from  that  assumed  by  Bekker  because  a  roil  particle 
moving  along  the  wheel  has  velocity  gradients  depending  again 
on  the  whole  flow-field  and  it  starts  to  be  sheared  before  I o 
reaches  the  wheel  surface. 

Concluding  these  remarks  on  Bekker' s  theory  it  may  be  said 
that  in  spi^e  of  the  progress  marked  by  this  theory  in  the  quan¬ 
titative  representation  of  the  soil-wheel  interaction,  it  has 
not  clarified  the  mechanism  of  soil  flow  beneath  the  wheel. 

Uffelmann  (l8)  has  assumed  that  the  normal  stress  on  a 
wheel  is  constant  and  equal  to  that  under  a  static  rough  foot¬ 
ing  in  a  plastic  Incompressible  clay  and  he  suggested  the  com¬ 
putation  of  the  wheel-resistance  on  this  basis.  Again,  because 
of  the  different  plastic  flow  patterns  in  the  two  cases,  the 
pressures  are  generally  different.  The  normal  stress  is  con¬ 
stant  on  a  footing,  but  varies  along  a  wheel. 

Reece  and  Wong  (13)  have  refined  Bekker 's  theory 
by  empirically  relating  the  position  of  the  point  of  maximum 
normal  stress  on  the  wheel  to  slippage.  Consequently  the  well- 
known  fact  that  the  tangential  stress  influence  tne  distribu¬ 
tion  of  normal  stresses  is  reflected  by  the  formulae  (fact 
Ignored  by  both  Bekker  and  Uffelmann  theories).  Again,  this 
generalization  is  based  on  some  empirical  data  and  not  on  the 
analysis  of  the  flow-field. 
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In  conclusion.,  there  is  no  comprehensive  theory  of  soil- 
flow  beneath  a  wheel.  All  existing  approaches  are  empirical  and 
consider  the  soil-wheel  interaction  as  a  surface  effect.  They 
ignore  the  stress  and  velocity  distribution  in  the  soil  body.  No 
consideration  has  been  given  to  inertial  or  strain-rate  effects. 

B.  Experimental  Work 

The  experimental  works  on  soil-wheel  interactions  are  dom¬ 
inated  by  the  same  trend  as  the  theoretical  works:  only  forces 
acting  on  wheels  or  stress  distribution  have  been  measured. 

Experiments  in  which  the  flow-pattern  of  the  soil  has  been 
observed  are  reported  in  (21),  (22),  most  of  them  being  carried 
out  with  sand.  Because  of  the  small  dimensions  of  the  experi¬ 
mental  soil  bin  and  the  low  speed  of  the  wheel,  the  quantitative 
value  of  the  results  is  questionable.  Moreover,  no  stress  measure 
ments  -  to  be  correlated  with  the  kinematics  -  have  been  carried 
out.  Qualitatively,  however,  these  experiments  confirm  beyond 
any  doubt  that  the  dynamics  of  wheel  motion  is  determined  by 
the  soil  flow  beneath  the  wheel  and  not  Just  by  a  surface  effect 
at  the  interface. 

The  stress  distribution  along  wheels  and  the  resulting 
forces  have  been  measured  by  different  authors,  but  there  is  no 
general  agreement  between  measurements .  For  instance  Uffelmann 
(l8)  has  found  that  the  maximum  normal  stress  acts  at  the  bottom 
of  the  wheel,  while  in  (15)  and  (13)  it  has  been  found  somewhere 
forward  near  the  bow.  Any  analysis  is  hampered  by  the  fact  that 
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results  are  dependent  on  the  type  of  soil  and  Its  preparation. 
Since  generally  there  is  no  detailed  description  of  the  mechan¬ 
ical  behavior  of  the  soil  used  by  different  experimenters 
it  is  difficult  to  compare  different  results  on  a  common  basis. 

The  comprehensive  measurement  of  forces  on  towed  and 
driven  wneels  by  Cullen,  Cullingford  and  Mayfield  (5)  gives  some 
qualitative  trends  as  the  decrease  in  flotation  with  slippage 
at  constant  sinkage. 

Some  experiments  are  compared  with  the  results  of  the 
present  work  in  Section  9* 

In  conclusion,  there  is  an  urgent  need  for  careful  experi¬ 
ments  in  which  both  soil  velocity  distributions  and  stress  dis¬ 
tributions  along  the  wheel  should  be  measured.  These  experiments 
have  to  be  carried  out  under  thorough  control  of  the  soil  be¬ 
havior  and  respecting  the  inertial  similitude  requirements. 

C ,  Theoretical  Works  or;  Rigid-Plastic  Flow 

Although  the  present  study  relies  on  the  rigid-plastic 
theory,  it  is  beyond  our  purpose  to  review  the  publications  con¬ 
sulted  for  this  study.  The  most  referred  work  is  Hill's  book 
(8)  and  to  a  lesser  extent  Prager's  (10),  (ll)  and  Thomas  (17) 
books  The  results  of  some  papers  of  special  interest  (l),  (9)> 
are  recalled  in  some  detail  in  the  context  of  the  following 
sect  1 ons . 
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3.  THE  STATEMENT  OF  THE  PROBLEM  AND  THE  GENERAL 
EQUATIONS  OF  RIGID- PLASTIC  FLOW 

A .  The  Statement  of  the  Problem 

Under  the  assumptions  enumerated  in  Section  1,  the  soil- 
wheel  problem  may  be  stated  in  the  following  terms  (Figure  l): 

A  standing  rigid  cylinder  of  radius  r'  rotates  with  con¬ 
stant  speed  u) '  and  has  a  given  constant  sinkage  z'.  The  soil 
approaches  the  wheel  with  a  constant  speed  at  infinity  U  '  and 
unperturbed  free-surface  elevation.  Find  the  free-surface 
shape,  the  soil  velocity  and  stress  fields  in  the  whole  semi- 
infinite  domain  and  particularly  along  the  wheel.  The  normal 
and  shear  stresses  along  the  wheel  (o'  and  T')  result  in  a 
vertical  force  W1 ,  a  horizontal  force  H1  and  a  torque  M1 . 

The  equations  of  rigid  plastic  flow,  on  which  the  solu¬ 
tion  is  based  in  this  study,  are  discussed  in  the  following 
paragraphs.  The  boundary  conditions  are  analyzed  in  Section  5* 

The  approximate  solution  is  given  in  Sections  6  and  7. 

In  applications  the  sinkage  z'  is  not  given,  but  the  verti¬ 
cal  load  W' .  From  a  mathematical  point  of  view,  however,  it  is 
much  more  convenient  to  start  with  a  given  sinkage.  Once  the 
problem  is  solved,  a  correlation  between  sinkage  and  flotation 
is  found  and  z'  may  be  determined  for  a  given  W' . 

B.  The  Physical  Behavior  of  a  Rigid-Plastic  Material 

The  ideal  rigid-plastic  model,  adopted  here  for  representing 
the  clay  behavior,  is  based  on  the  following  assumptions  (8),  ( 17 ) 
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-  The  plastic  deformations  are  much  larger  than  the  elas¬ 
tic  ones  (which  may  be  neglected).  The  elastic  state  is  there¬ 
fore  approximated  by  a  rigid  state  (Young  modulus  tends  to  in¬ 
finity).  The  only  possible  rigid  motions  are  translation  and 
rigid  rotation. 

-  In  the  plastic  state  the  maximum  shearing  stress  at  any 
point  is  constant  and  equal  to  the  yj.eld  stress  k'  .  In  the 
rigid  state  T '  <  k '  . 

-  The  material  in  plastic  state  flows,  the  direction  of 
maximum  shear  stress  and  the  direction  of  maximum  shearing 
strain-rate  being  parallel  (i.e.  the  stress  and  strain-rate  ten¬ 
sors  have  parallel  principal  axes).  There  is  no  one  to  one 
correspondence  between  the  stress  and  the  strain-rate,  which 
means  that  in  a  plastic  flow  the  maximum  shear  stress  is  con¬ 
stant  but  the  shearing  strain  rate  may  vary.  The  deformation 
work  is  totally  irreversible  in  a  thermodynamic  sense. 

-  The  material  is  incompressible  and  a  hydrostatic  pres¬ 
sure  does  not  influence  the  plastic  flow.  Only  the  devlatorlc 
stress  plays  a  role  in  creating  plastic  flow.  The  material  is 
homogeneous  and  isotropic. 

The  above  assumptions  lead  to  the  simplest  rigid-plastic 
model.  They  may  be  modified  in  order  to  take  into  account  elas¬ 
tic  effects,  anisotropy,  friction  and  work-hardening.  Although 
this  study  is  restricted  to  the  application  of  the  ideal  rigid- 
plastic  model,  it;  is  worthwhile  to  discuss  briefly  the  work- 
nardening  model,  which  is  useful  for  future  work  on  wheel -soil 
interaction. 
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While  an  ideal  non-hardening  material  yields  under  constant 
yield  stress,  otherwise  being  rigid,  the  yield  stress  of  work¬ 
hardening  material  is  dependent  on  the  history  of  strain  and 
increases  with  the  total  shearing  strain  (or  with  the  work  of 
deformation).  The  stress-strain  curve  of  a  non-hardening  ideal 
rigid-plastic  material  is  a  straight  line  (Figure  2),  while  in  a 
work-hardening  material  it  is  curved.  The  physical  microscopical 
mechanisms  which  underlie  work-hardening  of  metals  (8)  and  of 
soils  are  probably  different,  but  the  macroscopical  description 
is  satisfactory  in  both  cases.  Due  to  mathematical  difficulties 
the  work-hardening  theory  has  apparently  not  yet  been  used  in 
soil  mechanics.  An  ideal  rigid-plastic  material  (meaning  also 
non-hardening  here)  may  be  regarded  as  a  limit  of  a  work-harden¬ 
ing  material  with  a  very  high  rate  of  hardening  (Figure  2).  This 
interpretation  will  be  useful  in  the  discussion  of  Section  4c. 

C .  The  Equations  of  Two-Dimensional  Steady  Flow  (Cartesian 
Coordinates ) 

Although  the  equations  of  plastic  flow  may  be  found  in 
text-books  (3),  it  is  worthwhile  to  discuss  them  here  briefly. 


The  kinematics  of  flow  is  described  with  the  aid  of  the 


velocity  components  u',  v',  and  the  stress  field  by  a  ',  a  1 

“  y 

Txy'  (Figure  3). 


are 


The  equations  of  motion,  valid  for  any  type  of  material 
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a  c 1 
x 

ax' 
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dT  ' 

_ EL  = 

ay' 


p' 


au' 
3x  1 


+  V1 


du  1 

ay' 


[3.1] 


a  t  '  ao' 

_ EL  +  _ L  = 

ax'  ay1 


p' 


Sv ' 

T~T  +  v 

a>: ' 


ay 1 

ay' 


+  P'g 


C3.2] 


The  posit ive  signs  of  the  different  quantities  being  those  of 
Figure  3 • 

The  assumption  of  incompressibility  yields  the  continuity 
equation 


3u  1  ay 1 
ax'  ay' 


C3.3] 


The  assumption  of  constancy  of  yield  stress  in  the  plastic 
zone  is  reflected  by  Von-Mises  criterion  (8) 


(o  '  -  a  '  f  +  4  t  '  3  =  4k'3 

v  x  y  xy 


[3.4] 


and  finally,  the  parallelism  of  principal  directions  of  stress 
and  strain-rate  tensors  is  mathematically  described  by  Saint- 
Venant  relationships 
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which  by  the  continuity  equation  become 


a  1  - 
x 


2  T  ' 


■  a 
__Z. 

xy 


i 


3u  '  5 y  1 

5x'  5y 1 
Su1  av^ 
ay  •  ax' 


[3.6] 


The  five  differential  Equations  [3.l]>  [3.2],  [ 3 . 3 ] »  [3.^]  and 

[3.6]  are  in  principle  sufficient  in  order  to  determinate  the 

unknown  functions  a  o  T'  ,  u'  and  v!. 

x  y  xy 

The  five  equations  may  be  reduced  in  a  particularly  con¬ 
venient  way  to  a  system  of  four  equations  by  introducing  Von  Mis 

variables  p'  =  -  (a  1  +  a  ')/2  (the  isotropic  pressure)  and  0-the 

x  y 

angle  between  the  direction  of  maximum  shearing  stress  and  the 
x'  axis.  The  two  orthogonal  directions  of  maximum  shearing 
stress  a  and  £  (Figure  4)  are  selected  such  that  the  maximum  al¬ 
gebraic  normal  stress  acts  in  the  first  or  third  quadrant  (Fig¬ 
ure  4). 

The  new  variables  p',  0  replace  ' ,  a  '  and  through 
the  relationships 

a  '  =  -p'  -  k'  sin  2  0 

X 

a  1  =  -p1  +  k'  sin  2  0  [3.7] 

y 

T 1  =  k'  cos  2  0 

xy 


which  satisfy  Von  Mlses  criterion  [3.^]. 
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Using  dimensionless  variables  (defined  in  "Notation")  the 
equations  of  plastic  flow  become 


^  +  2  cos  20  ~  +  2  sin  20  ~  = 


dx 


ay 


-  H 


^  +  2  sin  20  -  2  cos  20  —■  =  -  H 


dx 


ay 


C3.91 


cos  20 


du  dv  ,  ,  du  ,  dv 

r—  -  r—  +  sin  20  r —  +  r — 
OX  oy  dy  dx 


=  0 


C3.10] 


in.  +  =  0 

dx  dy 


[3.11] 


1  II  |3 


The  plastic-inertial  dimensionless  number  H  = 


p'U 


(some¬ 


how  similar  to  Reynolds  number  in  a  viscous  flow)  is  particularly 
significant  for  the  soil-wheel  interaction. 


D.  The  Characteristic  Directions  and  the  Inertial  Terms 


In  plasticity  theory  the  inertial  and  gravitational  terms 
are  neglected  (equivalent  to  assuming  H  =  O.G  =  0).  In  this 
case  Equations  [3-8]  -  [3. 11]  form  an  uncoupled  quasi-linear 
system  of  first  order:  the  derivatives  of  p  and  0  appear  in  the 
two  first  equations,  while  those  of  u  and  v  in  the  last  equations. 
Each  of  the  two  systems  is  totally  hyperbolic  and  has  two  charac¬ 
teristic  directions  (8).  Accidentally,  the  characteristic  direc¬ 
tions  for  p,0  and  u,v  coincide  and  are  the  a  and  j3  directions. 
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The  equations  are  rewritten  in  curvilinear  characteristic  co¬ 
ordinates  x  ,  xQ  by  using  the  variables  transformation  (8) 

ct  p 


a 

ax 


_a_ 

dy 


V 

a 


cos  0 


sin  0 


a 

u  cos  0  + 
-  u  sin  0 


+ 


v 


+ 


sin  0 

cos  0 

sin  0 

v  cos 


a 

a 


0 


[3.12] 


the  velocity  being  now  referred  to  the  new  coordinates. 

After  some  algebraic  manipulations  the  system  [3.8]  -  [3.1l] 
becomes  (with  H  =  G  =  0) 


dx 
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dp  p  _££. 
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[3.13] 
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C  3 . 1^3 
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Equations  [3. 1^3  are  known  as  Geiringer  equations.  The 
lines  a  =  const  and  £  =  const  are  known  as  sliplines  and  they 
determine  an  orthogonal  network. 

Equations  [3.133  and  [3.143,  known  as  the  equations  of  quasi¬ 
static  plasticity,  have  the  well-known  properties  of  hyperbolic 
systems:  The  Cauchy  problem  for  p,  0,  v  ,  vQ  has  a  unique  solu- 

ct  p 

tlon  if  the  data  are  given  on  a  noncharacteristic  curve;  the 
solution  at  a  point  depends  only  on  the  data  on  the  portion  cut 
by  the  two-characteristics  through  the  given  point;  velocity 
tangential  discontinuities  are  possible  along  n  characteristic; 
the  velocity  derivatives  or  p  and  0  derivatives  may  be  dlncontin- 
uous  along  charac  eristics.  p  and  0  may  be  dl scontl nuouu  n  1  ong 
a  line  of  stress-discontinuity  which  cannot  br  i  chnrn  terlatl 
All  the  above  properties  are  discussed  in  detail  by  HI  , 

and  will  be  used  extensively  in  Sections  5  and  f  , 

The  neglectlon  of  the  inertial  terms  In  p  n  >M  U,  t 
when  applied  to  metal  technology  i3  entirely  Ji  I ' I 
H  «  1  there  (because  of  the  very  high  vnlu  I  , 
is  true  for  foundation  engineering  whorr  t 
equilibrium  are  considered.  In  the  cane  !  wb 
tion,  however,  H  may  be  of  the  order  0(1).  It 

a  soil  with  k'  =  3  psl,  p'g  =  120  pcf,  H  L» 
for  a  trans lational  speed  U^'  =7.5  mph.  h<  It  I  1 
therefore,  may  be  important  in  mobility  pr  >1  i>  m 
has  been  ignored  in  both  plasticity  theory  nr,  1  1  !  1  > 

Most  of  the  laboratory  tests  have  been  cnrrl  1 
of  H  much  smaller  than  in  field  conditions. 
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An  attempt  to  discuss  the  full  equations  of  plastic  flow  is 
presented  by  Spencer  (l6).  Spencer  analyzes  the  system  C 3 . 123 
as  a  quasi-linear  system  for  p,  0,  u,  v  and  determines  the  char¬ 
acteristic  directions  by  the  well-known  method  (see  (4),  p.  171- 
173). 


The  characteristic  determinant  of  the  system  [3.8]  -  C 3 . Ill 
considered  by  Spencer  has  the  form 


1  cos  20  - 

-  v  sin  20 

-X  sin  20  ♦ 

*  \  c  oe  20 


H(u  -  Xv) 


cos  20  - 
-  v  gin  20 


H(u  -  Xv) 


sin  20  + 

+  X  cos  20 

-X 


=  0 


1  h  i  .  !  111  11  reel  Ion . 

'  *  il  /i M  h  w,  that  rather  than 

•  i  |  y  i w  r  tip  r  X 


t  tan  0  [3.15] 

. »  *  H 

■  4-4#  f*t* 

t  ►  I  1 1  11 


mu  lorlatl  of  the  quasl- 
1  1  fir  that  the  system  [3.8]  - 

1  arc  11  ha  mcleristics .  This 
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conclusion  is  open  to  criticism  for  the  following  reason:  Like 
in  the  quasi-statical  case,  Equations  C 3 . 103  and  [3.1l3  form  a 
decoupled  system  for  u  and  v  which  is  hyperbolic  and  has  the 
characteristic  directions  of  [3*153.  No  matter  what  the  form 
of  the  right  hand  sides  of  Eqations  [3.8]  and  [3*93,  the  deter¬ 
minant  will  vanish  for  these  two  characteristic  values  of  X, 
which  can  be  easily  seen  from  the  inspection  of  the  last  two 
lines  of  the  determinant.  No  conclusion  regarding  the  Influence 
of  the  right  hand  side  of  Equations  [3.83  and  [3*93  can  be  drawn 
from  its  vanishing.  A  better  insight  into  the  problem  may  be  ob¬ 
tained  by  writing  the  system  [3.83  -  [3.113  in  characteristic 
coordinates  i.e.,  by  the  aid  of  Equations  [3.123.  The  result  is 
(16) 


dp  0  30  ( avq  30  ' 

3™—  +  2  r -  =  -  H  V_  r -  -  V„  r - 

SXa  M  P  ^  H 


*£._  .  Hv  r;i+v  M 

3 x„  3xn  a  3x  a  3x 

j3  /3  a  a 


+  G  sin  0 


>  [3.163 


+  G  cos  0 


3v 

a  30 

3x  “  VjS  3x 
a  r  a 


Svj3  30 

4.  v  - - 

Sxp  a  3xp 


=  0 


[3.173 


=  0 
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Obvlously  the  Geirlnger  Equations  C 3  - 17 3  are  the  same  as 
[3.143.  Equations  [3.16],  however,  do  not  satisfy  the  require¬ 
ment  of  fully  hyperbolicity  ((4),  p.  173)  which  states  that  all 
the  quantities  should  be  differentiated  along  the  same  charac¬ 
teristic.  The  presence  of  the  cross-derivatives  in  the  right- 
hand  sides  of  Equations  C3.I6]  shows  that  the  usual  properties 
of  a  hyperbolic  system  are  not  met  here. 

It  is  beyond  the  purpose  of  this  study  to  investigate  the 
nature  of  the  system  [3.163  and  [3.173  which  in  fact  is  the 
basis  of  a  new  chapter  of  applied  mechanics  which  can  be  called 
dynamic  plasticity.  At  this  stage,  assuming  that  for  the  soil- 
wheel  interaction  H  <  1,  a  perturbation  expansion,  similar  to 
that  suggested  by  Spencer  (l6),  will  be  adopted.  In  such  an 
expansion,  in  which  H  is  a  small  parameter,  the  leading  term 
satisfies  the  hyperbolic  quasi-static  equations  [3*133  and 
[3.143.  The  wheel-problem  will  be  studied  (Sections  5>6)  in 
this  first  quasi-static  approximation  at  this  stage.  In  Sec¬ 
tion  4  the  Influence  of  the  inertial  terms  on  singular  points 
and  lines  will  be  discussed  and  the  quasi-static  solution  will 
be  built  in  the  light  of  the  results.  In  this  sense  the  influ¬ 
ence  of  the  inertial  terms  will  be  reflected  in  the  proposed 
solution,  which  is  different  from  the  type  of  solutions  adopted 
in  rolling  and  extrusion  theories  (8),  (l). 

It  is  worthwhile  to  mention  here  that  in  a  hypothetical 
case  of  very  high  H  (for  Instance,  for  a  vehicle  which  is  not 
wheel-propelled)  a  perturbation  expansion  of  Equations  L 3 . l63 , 
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[3.17]  starts  with  the  equations  of  inviscid  flow  (Euler  equa¬ 
tions)  as  the  leading  term.  The  soil  behaves  in  this  case  as 
an  Inviscid  fluid,  rather  than  a  plastic  solid.  The  Euler  equa¬ 
tions  have  an  elliptic  character.  The  problem  for  Intermediate 
H  is  therefore  somewhere  between  these  two  extremes. 

In  an  analysis  of  the  equations  of  flow,'  written  in  the 
form  of  a  second  o"der  system  for  p,  u,  v  by  the  elimination  of 
0  from  13-83  -  C 3 •  111 ,  Thomas  (17)  has  found  the  same  character¬ 
istic  directions  L  3  - 1 5 3 ,  but  he  did  not  discuss  at  all  the  nature 
of  Equations  [3.16],  [3.171.  Assuming  continuous  velocities  and 
pressures,  Thomas  showed  that  discontinuities  in  the  second 
derivatives  of  the  velocities  can  occur  only  along  the  charac¬ 
teristics.  The  only  applications  discussed  by  Thomas  are  of  un¬ 
steady  flows  with  no  inertial  convective  terms. 

E.  The  Equations  Along  Streamlines 

The  equations  of  plastic  flow  have  been  written  so  far  in 
cartesian  and  characteristic  coordinates.  It  is  worthwhile  to 
rewrite  them  in  coordinates  attached  to  the  streamlines. 

With  0  the  angle  between  the  streamline  and  x  axis,  s  and 
n  coordinates  along  and  normal  to  the  streamline  (Figure  3)>  the 
following  relationships  hold 


=  _L 

9x  9  s 

_9_  _  _9_ 
9y  9  s 


cos  9 

sin  0 


-  r—  sin  0 
9n 

+  ^  o°s  9 


u  =  V  cos  0 
v  =  V  sin  9 


[3.18] 
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The  substitution  of  Equation  [3.18]  into  Equations  [3*8]  - 
[3.11]  gives 

|^+2  cos  2(0-0  )  |^  +  2  sin  2(0-0  )  =  -  H  -  G  sin  0 

^  +  2  sin  2(0-0  )  -  2  cos  2(0-0  )|^=-HVa  -  G  cos  0 

cos  2(0-0  )  ^  +  sin  2(0-0  )  ~  +  V  sin  2 (0-0  )  - 

-  V  cos  2(0-0  )  P-  =  0 
on 


51  +  v  —  - 

3s  3n 


[3.193 


The  Equations  [3. 193  are  used  in  Section  4  for  some  applications. 

F .  The  Condition  of  Positive  Work 

The  Saint-Venant  condition  [3.63  expresses  the  parallelism 
between  the  principal  directions  of  stress  and  strain-rate  ten¬ 
sors,  but  does  not  specify  the  signs  of  the  two,  which  permits 
for  two  alternatives.  Only  one  of  them,  which  satisfies  the  re¬ 
quirement  of  positive  dissipation  work  is  a  valid  representation 
of  a  plastic  flow. 

The  power  of  dissipation  is  given  by 


E  =  e  a  +  2T  y 

xx  x  xy 'xy 


+  e  a 


yy  y 


>  0 


[3 .203 
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In  terms  of  u,v  and  0,E  is  expressed  as 


E  =  sin  20 


\ 


3v  3u  ,  du  ,  Sv 

r—  -  r—  f  cos  20  r-  +  7- 
oy  ox  dy  Ox 


>  0 


[3.21] 


In  characteristic  coordinates,  Equation  C 3 . 2 1 3  transforms 
with  the  aid  of  Equations  [3. 12]  into 


v  v0 
a  jB 


0V‘ 

dx 


+  v 


3V‘ 


a 


/3  3x 


0 


>  0 


[3.22] 


Equation  [3. 22]  has  been  found  by  Prager  (10)  from  graphical  con- 

d  d 

siderations.  Since  the  operator  vq  +  v_  represents  dif- 

a  xa  P 
.  a . 

ferentlatlon  along  streamlines  (7—),  E  may  be  finally  expressed 

0  s 

as 


E  = 


v  v_ 
a  0 


>  0  . 


[3.23] 
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4.  SOME  APPLICATIONS  OF  THE  GENERAL 
EQUATIONS  OF  PLASTIC -RIGID  FLOW 

As  stated  in  Section  3>  the  general  study  of  the  equations 
of  plastic-rigid  flow  [3.8]  -  [3.  Ill  or  [ 3 . 16] ,  C3.17]  is  be¬ 
yond  the  scope  of  this  report.  A  discussion  of  the  application 
of  the  equations  to  some  particular  lines  and  singular  points 
is  presented  here  because  of  its  importance  for  the  wheel  problem. 

A .  The  Free-Surfaoe 

The  soil  free-surface  in  the  rigid  zone  of  uniform  flow  is 
horizontal.  In  the  vicinity  of  the  wheel, in  the  plastic  region, 
the  free-surface  curves.  Three  boundary  conditions  have  to  be 
satisfied  on  a  free-surface 


a  =  0  >> 

} 

[4.1] 

T  =  0  J 

n 

V  =  0 
n 

[4.23 

i.e.,  the  free-surface  is  stress-free  and  it  is  a  streamline  (in 
a  steady  flow). 

From  Equations  [3. 73  we  immediately  find  that  Equations  [4.l] 
are  equivalent  to  (Figure  5) 

0  -  0  -*  J 


p  =  ±  1 


[4.33 
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i.e.,  the  a  lines  intersect  the  free-surface  at  45  degrees.  The 
plus  sign  in  Equations  C  4 . 3 1  corresponds  to  a  state  of  compres¬ 
sion,  the  maximum  compressive  stress  being  a  =  -2,  For  0  -  0  =  -tt/4 
the  normal  stresses  are  tensile  with  a  maximum  value  of  j  =  2,  i.e. 
twice  the  yield  stress  (Figure  5). 

Considering  first  the  state  of  compression,  it  is  easy  to 
ascertain  that  (Figure  5a) 


v 

a 


v 

V  cos  t- 


v 


0 


-V  sin  £ 


[4.4] 


Hence,  from  the  condition  of  positive  work  [3.231  it  is  found 
that 


3Va 

3s 


< 


0 


[4.51 


which  means  that  the  speed  along  a  free-surface  has  to  decrease 
in  the  direction  of  flow. 


The  equations  of  flow  [3-191  become  in  the  case  of  a  free- 
surface  (0  -  0  =  tt/4 , p  =  l) 

„  90 


3n 


H  av^_ 

2  3s 


-  G  sin  0 


+  2  —  = 
3n  3  s 


-  HV3  P-  -  G  cos  0 
0  s 


dn  3s 


=  0 


T—  +  V  r—  = 
3  s  3n 


[4.6] 


0 
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The  second  equation  of  [4.6]  shows  that  for  a  convex  free-surface 
the  pressure  decreases  along  the  normal,  inside  the  flow  region, 
i.e.,  tensile  stresses  are  Induced  beneath  the  free-surface.  The 
last  equation  of  [4.6]  shows  that  90 /9n  >  0,  i.e.,  the  stream¬ 
lines  diverge  in  the  flow  direction. 

In  the  case  of  tension  (0  -  0  =  -  tt/4,  p  =  l)  the  speed  has 
to  Increase  along  the  free-surface. 

B.  A  Boundary  With  a  Rigid  Body 

In  plasticity  theory  (8)  two  conditions  are  imposed  on  a 
boundary  with  a  rigid  body:  the  normal  speed  is  zero  (streamline) 
and  a  stress  condition.  Three  types  of  stress  conditions  are 
assumed:  on  a  "rough"  surface  T  =  1,  on  a  "smooth"  lubricated 

surface  T  =  0  and  on  a  frictional  surface  t  is  proportional  to 
o  (Coulomblan  friction).  There  is  no  restriction  on  the  tan¬ 
gential  velocity  which  is  generally  different  from  that  of  the  rigid 
body  (Figure  6).  The  sign  of  the  shearing  stress  on  the  rigid 
body  depends,  however,  on  the  relative  flow  along  the  body:  t 
acts  in  the  direction  of  relative  flow. 

Assuming  for  the  time  being  (see  Section  5)  that  the  stress 

boundary  condition  is  T  =  T  ,  i.e.,  the  shearing-stress  is  given 

w 

(0  £  t  £  l),  let  us  first  consider  the  Equations  [3.7].  For  a 
w 

given  t  the  angle  between  the  a  line  and  the  body  is  given  by 
w 

(Figure  6) 


T  =  cos  2(0  -  0 ) 
w 


[4.7] 
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Since  T  acts  in  the  direction  of  relative  flow,  the  angle 
w 

between  the  a  line  and  the  boundary  is  smaller  than  tt/4  with 
respect  to  the  relative  flow  direction,,  Three  cases  of  a 
locally  slipping  and  skidding  boundary  are  represented  in  Fig¬ 
ure  6. 

The  velocity  components  on  the  characteristic  directions 

are 

v  =  V  cos  (0  -  0 ) 

a  [4.8] 

v^  =  -V  sin  (0  -  8  ) 

The  condition  of  positive  work  C 3  -  23 3  applied  along  the 
boundary  gives 


1 

cos  (0  -  0)  sin  (0  -  9) 


ay3 

ds 


>  0 


[4.91 


Hence,  for  an  accelerating  flow  along  the  boundary  in  the  direc¬ 
tion  of  relative  flow  Equation  [4.91  gives 

sin  (0  -  0 )  >  0  ,  [4.10] 

since  cos  (0  -  0  )  >  0  for  1 0  -  0  |  *  jj-  .  Consequently,  in  the 
case  of  an  accelerating  flow  0  >  0  (Figure  6a),  while  for  a 
retarded  flow  0  <  0  (Figure  6b). 
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A  limit  case  is  particularly  important  since  it  permits 

the  exact  integration  of  the  equations  of  flow,  namely  T  =  1 

w 

(a  rough  boundary)  i.e.,  0  -  '9  =0. 

In  this  case  the  boundary  is  a  characteristic  line  and 
Equations  [3.191  become 


3p  30 

3  s  3  s 


H  3Va 

2  is 


G  sin  0 


p  90 

3n  3n 


H  V3  P-  -  G  cos  0 
3s 


.  v  o 

3s  3n 


3V 

3s 


+ 


0 


[4.11] 


From  the  last  two  equations  we  find  =  0,  V  =  const. 

Hence,  the  speed  is  constant  along  a  streamline  which  is  a  char¬ 
acteristic  . 

The  first  equation  of  [4.1l]  yields 

p  +  20  -f  Gy  =  const.  [4.12] 

Equation  [4.12]  shows  how  the  pressure  variation  depends  on 
the  curvature  of  the  boundary.  In  the  case  of  a  wheel  p  decreases 
along  the  boundary,  since  30 /3s  >  0  (assuming  that  the  gravity 
effect  is  negligible.)- 
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In  the  case  in  which  a  streamline  is  an  envelope  of  char¬ 
acteristics  the  velocity  may  vary  along  it.  This  property  has 
been  used  in  Prandtl  solution  (8)  of  flow  between  two  rigid 
plates  and  Alexander  (l)  solution  of  rolling.  This  possibility 
is  used  in  this  work  too  (see  Section  6). 

C .  A  Line  of  Separation  Between  Regions  of  Rigid  and 
Plastic  Flows 

A  body  moving  in  a  rigid-plastic  material  generates  a 
region  of  plastic  nonuniform  flow  in  its  vicinity.  Since  at  in¬ 
finity  the  flow  is  uniform  and  in  a  rigid  state, there  is  a  transi¬ 
tion  between  the  two  regions.  In  all  problems  of  metal  technol¬ 
ogy  similar  to  the  soil-wheel  problem  (extrusion,  rolling,  wedge 
Indentation  (l),  (8),  and  solved  with  the  aid  of  the  quasi-state 
equations ), it  is  assumed  that  this  transition  is  realized  through 
a  line  of  tangential  velocity  discontinuity. 

Such  a  line  is  regarded  as  a  limit  of  a  narrow  zone  of 
plastic  flow  in  which  the  tangential  velocity  varies  rapidly 
(Figure  7),  (8).  In  this  layer  the  direction  of  maximum  shear¬ 
ing  strain-rate  is  parallel  to  the  layer  axis.  It  is,  therefore, 
assumed  that  the  discontinuity  line  is  a  characteristic  line. 

Along  such  a  line  the  tangential  velocity  jump  is  constant. 
This  is  simply  found  from  Geirlnger's  Equations  [ 3.171  (Figure 
7b).  Assuming  that  the  discontinuity  line  is  an  a  line,  from 
continuity  normal  to  the  line 

v  '  =  -U  '  sin  0 
P  o 


[4.13] 
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The  first  of  Equations  [3.171  gives 

r — =  -U'sin  0  z~~—  or  [4.14] 

ox  1  o  ox  ' 

a  a 

v  1  -  U  1  cos  0  =  const  [4.151 

a  o 

Since  the  component  of  the  rigid  uniform  flow  along  the  line 
is  U  '  cos  0,  Equation  [4.151  states  that  the  tangential  jump  is 
constant . 

r 

In  a  rigid -plastic  material  the  shear  stress  in  the  plastic 
narrow  zone  is  constant  and  equal  to  the  yield-stress,  T'  =  k'. 

Since  there  is  no  variation  in  the  shear  stress  across  a  thin 
layer  when  inertial  effects  are  neglected,  the  shear  stress  at 
the  boundary  between  the  plastic  layer  and  the  rigid  zone  is 
T'  =  k'  and  decreases  inside  the  rigid  zone. 

In  a  work-hardening  material  in  which  inertial  effects  are 
neglected  the  discontinuity  layer  has  to  diffuse  and  to  have  a 
finite  thickness,  since  the  shearing  stress  at  the  limit  of  the 
rigid  zone  must  be  smaller  than  the  ultimate  yield-stress.  If  ' 

the  zone  of  work-hardening  in  the  stress -strain  diagram  (Fig¬ 
ure  2b)  is  narrow,  the  transition  layer  may  be  thin  and  the  dis¬ 
continuity  line  is  still  a  valid  approximation. 

When  inertial  effects  are  taken  into  account  the  first 
question  to  be  asked  is  if  a  line  of  velocity  discontinuity  is 
possible.  Applying  the  momentum  equation  to  such  a  discontinuity 


[ 
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line  (Figure  7c)  it  is  found  that  a  shear-stress  discontinuity 
exists  across  the  line 


5T  ' 


T  ' 

plastic 


T  ' 

rigid 


p'U  1  sin  0  6V 
r  o 


or  in  dimensionless  variables 


[4. 16] 


6T  =  H  sin  0  5V  . 


[4.171 


If  5T  is  much  smaller  than  unity  (5T1  «  k')  or  0  =  0  (no  cross- 
flow),  inertial  Jump  may  be  neglected  and  the  conclusions  ob¬ 
tained  from  quasi-static  considerations  are  valid.  But  in  the 
case  of  soil-wheel  interaction  H  may  be  of  order  unity  and  the 
inertial  Jump  cannot  be  neglected.  The  important  conclusion  is 
that  (assuming  that  in  the  transition  layer  the  characteristics 
are  directed  along  the  layer),  the  inertial  effect  causes  the 
thickening  of  the  layer  in  an  ideal  rigid-plastic  material, 

since  otherwise  the  plastic  flow  will  start  with  T '  <  k' . 

r  max 

Hence,  a  velocity  discontinuity  is  not  acceptable  in  an  ideal 
rigid-plastic  material. 

If  work-hardening  is  taken  into  consideration,  then  a 
narrow  transition  zone  becomes  possible  again  when  inertial  ef¬ 
fects  are  considered,  since  plastic  flow  may  start  with  T'  <  k' . 
At  any  rate,’  since  T  '  >  0  (Figure  2b)  the  inertial  Jump  is  lim¬ 
ited  by  the  value 


BT  =  H  BY  sin  0  <  1  . 


L4.18] 
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In  a  similar  context  the  direction  of  the  velocity  discontinuity 
has  to  be  such  that 


6T  =  H6V  sin  0  2  0  [4.191 

i.e.,  the  maximum  shearing  stress  in  the  rigid  zone  is  smaller 
than  the  yield  stress. 

At  this  stage  it  will  be  assumed  that  in  the  case  of  soil- 
wheel  interaction  the  rigid  and  plastic  zones  may  be  separated 
by  a  line  of  velocity  discontinuity,  provided  that  Equations 
[4.18]  and  [4.191  are  satisfied.  It  is  however,  fair  to  say 
that  this  point  needs  further  theoretical  and  experimental  in¬ 
vestigation  in  the  framework  of  dynamic  plasticity. 

The  velocity  jump  is  constant  and  given  by  Equation  [4.15], 
since  inertial  terms  do  not  affect  the  Gelrlnger  equations. 

D .  The  Centered  Fan 

A  bundle  of  Intersecting  characteristics  of  the  same  family 
determine  a  centered  fan  (Figure  8a).  The  vertex  of  the  fan  is 
a  singular  point.  This  singularity  plays  an  important  role  in 
quasi-static  plasticity  and  it  is  used  in  the  solution  of  ex¬ 
trusion,  rolling  and  sheet-drawing  problems  (l),  (8).  In  all 
these  solutions  a  velocity  discontinuity  is  admitted  at  the 
vertex.  It  will  be  subsequently  shown  that  in  a  plastic  flow 
with  Inertial  effects  such  a  discontinuity  is  not  possible. 
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Assuming,  for  sake  of  definiteness,  that  the  straight  char¬ 
acteristics  are  a  lines,  we  may  replace  the  a,  /3  coordinates  by 
polar  coordinates  (Figure  8) 


d  a 


3  13 


3x  3r  *  9xq  r  96 
a  0 


!  va  =  vr  :  "  v9  C’+-19:1 


while,  by  definition 


S0_  _  S0  1 

3x  3xfl  r 

a.  0 


[4  20] 


Substituting  the  above  relationships  in  Equations  [3.16] 
[3.17]  we  obtain 


=  _ 
3r 


H 


r 


v„ 


G  sin  0 


30 


2 


H  v 


r 


3r 


G  cos  0 


[4.21] 

[4.22] 


3v 

r 

3r 


0 


30  +  Vr 


0 


[4.2.3] 


Equation  [4.21]  implies  that  generally  p  becomes  Infinite  at  tne 
vertex.  This  inadmissible  effect  may  be  avoided  in  three  cases: 
(i)  vQ  =  0  which  by  Equation  [4.23]  gives  vr  =  0,  i.e.. ,  rest; 
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3v 

(ii)  -rg-  -  Vq  =  0  which  together  with  Equation  [4.23]  provides 
for  vr  and  Vg  the  solutions 

vr  =  Ql  sin  0  +  Cs  cos  0  [4.24] 

Vg  =  Ci  cos  0  -  C3  sin  0 

In  a  cartesian  system  the  velocity  distribution  [4.24]  is  equiva¬ 
lent  to 


COS 

9  -  Vg 

sin 

0 

=  cs 

sin 

9  +  Vg 

cos 

0 

=  Ci 

[4.25] 


representing  a  uniform  flow. 

The  pressure  distribution  is,  according  to  Equations  [4.21] 
and  [4.22] 


p  =  20  -  Gy 


[4.26] 


which  is  obviously  identical  with  that  obtained  in  a  quasi-static 
case; 


(m) 


v9  =  vg (r) 


[4.27] 


■"•UTAH 
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Equations  [4.231  and  [4.26]  shows  that  v  =  0  and  this  is, 
therefore,  a  rotational  flow  around  the  vertex  The  pressure 
distribution  in  the  centered  fan  is  given  by 

f  vQa 

p  =  20  -  Gy  +  H  I  — —  dr  [4.283 

Case  (ill)  offers  the  possibility  of  matching  plastic  flows 
around  concave  corners  (Figure  8b).  In  the  case  of  a  convex 
corner*,  in  quasi-static  solutions  it  is  generally  assumed  that 
the  flow  passes  from  the  uniform  regime  to  the  plastic  flow  (Fig¬ 
ure  8c)  through  a  line  of  velocity  discontinuity  AD  followed  by 
a  centered  fan  with  discontinuous  ADC  velocity  at  the  vertex  A 
(l),  (8),  (10).  In  a  flow  with  inertial  terms  taken  into  ac¬ 
count,  the  only  possible  type  of  flow  behind  the  discontinuity 
line  is  a  uniform  flow  (Figure  8d).  This  important  difference 
has  to  be  kept  in  mind  when  solving  the  soil-wheel  problem. 

E.  A  Line  of  Stress  Discontinuity 

The  shear  (T)  and  normal  stresses  ( a )  have  to  be  con¬ 
tinuous  across  a  line.  The  normal  stress  parallel  to  the  line 
(a  )  may  be  discontinuous  (Figure  9)°  Hill  ((8),  p  157)  has 
analyzed  this  type  of  discontinuity  along  a  line  separating  two 
plastic  zones. 

His  conclusions  are:  (i)  such  a  line  cannot  be  a  character¬ 
istic;  (ii)  the  characteristics  of  the  same  family  from  the  two 
sides  intersect  the  line  at  a  same  angle  \  (Figure  9)  »  (Hi) 
the  pressure  Jump  is  given  by  |  pi  -  Pa  |  =  2  sin  A,  and  (iv)  the 
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velocity  is  continuous  across  the  stress  discontinuity  line.  The 
last  property  shows  that  such  a  line  is  acceptable  in  a  plastic 
flow  when  inertial  effects  are  taken  into  account  and  all  above 
conclusions  remain  valid. 
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5.  APPLICATION  TO  THE  WHEEL  PROBLEM 


A .  General 

A  rigid  cylindrical  wheel  (roller)  rotating  and  moving 
with  constant  speeds  is  considered  here.  The  study  is  limited 
here  to  a  driven  slipping  wheel  which  by  definition  has  peripheric 
rotational  speed  larger  than  the  soil  speed  at  any  point  of  the 
soil-wheel  interface.  Consequently  the  shearing  stress  along 
the  wheel  is  positive  everywhere  on  the  interface. 

Assuming  that  the  soil,  a  soft  undrained  clay,  behaves 
like  a  rigid-plastic  material, one  has  to  solve  mathematically 
the  problem  by  integrating  the  equations  of  plastic  flow  [3.8]  - 
[3.1l]  with  the  appropriate  boundary  conditions. 

Due  to  the  difficulties  caused  by  the  inertial  terms  a  per¬ 
turbation  expansion  technique  is  used  in  order  to  simplify  the 
equations  ('’6).  The  influence  of  the  inertial  terms  on  some 
particular  lines  (Section  4)  will  be,  however,  taken  into 
account . 

B.  The  Perturbation  Expansion 

Spencer  (l6)  has  suggested  a  perturbation  expansion  of  the 
equations  of  flow  C3.16],  [3.173  in  which  H  is  considered  a  small 
parameter.  Assuming  that  all  the  variables  may  be  expanded  in 
a  power  series  as 
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P  =  p0  + 

0  =  0  + 
o 


H  pi  +  Ha  pa  +  .  .  . 
H  0i  +  Ha0j  +  .  .  . 


+  H  v 

+  Hav  +  . . . 

oti 

Ofl 

+  H  vfl 

+  Hav.  +  . . . 

Pi 

C5.ll 


and  substituting  In  Equations  [3.16],  a  set  of  separate  systems 
of  equations  Is  obtained  for  each  order  of  approximation.  Here 
only  the  zero  order  terms  are  considered.  They  are  satisfying 
the  quasi-static  equations  C  3  - 13  3  and  [3.1^]>  assuming  that  G 
(gravity  effects)  is  also  small. 

Since  only  the  zero-order  approximation  will  be  subsequently 
considered,  the  zero  sub-index  Is  omitted. 

A  legitimate  question  is  that  of  the  validity  of  the  expan¬ 
sion  in  cases  in  which  H  is  not  small.  A  rigorous  investigation 
of  this  problem  is  suggested  for  a  later  stage  (Section  11). 
Intuitively,  it  seems  that  for  moderate  value  of  H  (H  <  l)  the 
inertial  effects  are  mainly  concentrated  in  the  discontinuity 
lines  rather  than  in  the  plastic  region.  It  is,  therefore,  hoped 
that  taking  certain  precautions  in  the  ze^o-order  approximation, 
the  zero-order  solution  is  an  acceptable  approximation  for  the 
soil-wheel  Interaction  at  moderate  speeds. 
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C .  General  Properties  of  Quasl-Statlc  Solutions 

The  solution  of  the  soil-wheel  problem  In  the  quasl-statlc 
approximation  Is  still  very  difficult,  as  may  be  ascertained 
from  the  examination  of  similar  solutions  In  metal  technology. 

In  fact,  almost  all  solutions  obtained  so  far  in  plasticity 
theories  are  approximate,  or  "incomplete."  A  complete  solution 
of  a  plasticity  problem  requires  the  determination  of  p,  0,  u,  v 
satisfying  the  equations  of  flow  [3.133  -  [3.1*0  in  the  plastic 
region  and  the  boundary  conditions,  and  the  determination  of 
a^,  and  in  the  rigid  zone  satisfying  the  elasticity  equa¬ 
tions  and  matching  the  stresses  at  the  rigid-plastic  interface. 
The  equations  are  elliptic  in  the  elastic  zone  and  hyperbolic 
in  the  plastic  zone.  The  location  of  the  line  of  separation  be¬ 
tween  the  two  zones  is  unknown  and  the  same  is  true  for  the 
free-surface,  if  it  occurs.  Since  the  complete  solution  is  very 
difficult  the  usual  approach  in  plasticity  is  to  look  for  ap¬ 
proximate  incomplete  solutions  which  give  lower  and  upper  bounds 
of  the  exact  solutions  (8),  (ll). 

A  lower  bound  is  found  by  assuming  a  plastic  pattern  which 
satisfies  only  the  stress  boundary  conditions,  while  kinematical 
boundary  conditions  are  neglected.  In  other  words,  If  a  stress 
distribution  In  both  the  assumed  plastic  and  rigid  zones  which 
equilibrates  the  external  forces  is  found,  then  the  exact  solu¬ 
tion  will  allow  for  larger  external  forces. 
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The  upper  bound  is  found  by  assuming  a  plastic  flow  pattern 
which  satisfies  the  velocity  boundary  conditions.  The  forces 
become  in  this  case  an  upper  bound  to  allowable  external  forces. 
The  solution  is  found  generally  by  assuming  a  priori  a  certain 
qualitative  distribution  of  discontinuity  lines,  centered  fans 
and  slip-line  field  and  trying  to  adjust  them  in  order  to  satisfy 
the  equations  of  flow  and  the  boundary  conditions. 

In  the  soil-wheel  problem  the  boundary  conditions,  for  a 
given  sinkage,  are  expressed  in  kinematical  terms  mainly  (see 
Section  5E).  It  is,  therefore,  natural  to  seek  an  incomplete 
upper  bound  solution.  Such  a  solution  may  be  eventually  proved 
to  be  an  exact  one  if  the  stress  field  may  be  continued  in  the 
rigid  zone.  Otherwise,  the  validity  of  the  solution  has  ulti¬ 
mately  to  be  checked  by  comparison  with  experiments  (see  Sec¬ 
tion  9). 

It  is  worthwhile  to  mention  here  that  in  view  of  the  dif¬ 
ferent  simplifying  assumptions  it  seems  somehow  futile  to  refine 
the  computations  beyond  a  certain  limit.  An  approximate  solu¬ 
tion  which  reflects  the  main  features  of  the  problem  and  shows 
the  characteristic  trends  should  be  satisfactory  in  the  soil- 
wheel  case. 

D.  Qualitative  Analysis  of  the  Problem 

The  solution  of  the  soil-wheel  problem  has  to  be  found  now 
by  assuming  a  slipline  field  which  has  to  satisfy  equations  of 
plastic  flow  and  boundary  conditions.  This  s  facilitated  by  a 
qualitative  analysis  of  the  problem. 
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First,  It  is  reasonable  to  assume  that,  similarly  to  all 
known  solutions  in  plasticity  literature,  the  plastic  region  is 
of  finite  extent  and  attached  to  the  body  (experiments  tend  to 
confirm  it,  at  least  at  low  speeds).  Moreover,  in  quasi-static 
conditions,  there  is  no  reason  to  assume  a  forward  influence  of 
the  driven  wheel.  This  is  the  case  at  least  in  the  solutions 
of  rolling,  extrusion  and  sheet-drawing  processes.  If  slight 
upstream  Influences  are  observed  in  experiments  (22),  it  may  be 
well  due  to  work-hardening  effects  (see  Section  9E)  not  taken 
into  account  in  a  rigid-plastic  model.  At  any  rate  in  all  tests 
of  driven  wheels  in  clays  no  sensible  upstream  disturbances  of 
the  free-surface  have  been  reported  (l8),  (5)  which  confirms  the 
assumption  of  upstream  rigid  flow. 

In  the  downstream  region  the  free-surface  has  to  leave  the 
wheel  smoothly  and  a  plastic  "wake"  of  finite  dimension  is 
created.  This  is  different  from  rolling  and  extrusion  processes 
in  which  it  is  assumed  that  the  free-surface  Is  horizontal  im¬ 
mediately  behind  the  body  (l),  (8),  by  locating  a  centered  fan 
and  a  discontinuity  line  at  the  rear.  Inertial  effects  do  not 
allow  for  such  a  type  of  solution  (see  Section  6c)  and  experi¬ 
ments,  even  at  low  speeds,  in  both  sand  and  clay  confirm  this 
conclusion  (18),  (21),  (22). 

Consequently,  the  free-surface  has  qualitatively  the  shape 
suggested  in  Figure  10. 

The  kinematics  of  flow  may  be  also  depicted  from  qualita¬ 
tive  considerations.  Since  the  flow  is  incompressible,  from 
continuity  requirements  the  streamlines  have  to  look  somehow 
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as  in  Figure  10a.  The  average  speed  in  the  plastic  zone  is  max¬ 
imum  at  the  bottom  of  the  wheel.  The  flow  is  converging  from  up¬ 
stream  to  the  bottom  section  and  diverges  downstream  in  order 
to  match  the  uniform-rigid  downstream  flow.  Obviously,  in  an 
incompressible  material  and  a  two-dimensional  flow  there  is  no 
rut  and  both  free-surface  and  speed  have  to  reach  the  unperturbed 
state  downstream. 

A  qualitative  picture  of  a  suggested  slipline  field  built 
in  the  light  of  the  above  observations  is  presented  in  Fig¬ 
ure  10b.  A  detailed  argumentation  and  examination  is  presented 
in  Section  £. 

At  the  bow  A  a  line  of  discontinuity  AB  separating  the  rigid 
and  plastic  zone  is  suggested.  This  is  an  a  line  since  shearing 
stresses  acting  on  the  rigid  material  have  to  have  the  same  sign 
as  those  acting  on  the  wheel.  The  speed  along  the  wheel  is  in¬ 
creasing  from  A  to  C,  so  that  the  j3  lines  have  to  intersect  AC 
at  the  angles  shown  in  Figure  10b  (according  to  Section  4b) .  On 
the  rear  portion  OR  it  is  assumed  that  the  flow  is  decelerated. 
Hence,  the  /3  sliplines  have  to  intersect  CR  at  the  angles  of 
Figure  10b.  The  transition  between  the  two  zones,  of  converging 
and  diverging  flow,  is  realized  through  a  rigid  core  BCD  which 
moves  horizontally  with  the  bottom  soil  speed.  A  somehow  simi¬ 
lar  rigid  region  is  assumed  in  rolling  solutions  (l).  Although 
not  necessarily  rigorous,  this  separation  of  the  two  plastic 
regions  ABCA  and  CBETRC  by  a  zone  of  uniform  flow  is  very  con¬ 
venient  for  computational  purposes  (Section  6).  The  zone  RTER 
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is  a  curvilinear  triangle  created  by  the  free-surface  TR.  The 
line  ET  separating  plastic  and  rigid  zones  at  the  rear  again 
has  to  be  a  characteristic  line  since  at  least  the  second  de¬ 
rivatives  of  the  velocity  are  discontinuous  there  (17).  This 
is  assumed  again  to  be  an  a  line,  with  the  rigid  zone  sheared 
in  the  same  direction  as  the  soil  beneath  the  wheel.  From  in¬ 
ertial  considerations  (Section  6C)  there  is  no  possible  veloc¬ 
ity  jump  on  the  line  ER  and  the  flow  has  to  match  smoothly  the 
uniform  downstream  flow. 

The  region  CDER  is  a  transition  region  in  which  the  flow 
diverges.  Along  DE  the  speed  drops  from  its  maximum  value  to 
the  unperturbed  speed.  The  only  possible  way  to  satisfy  this 
requirement  is  to  assume  that  DE  is  an  envelope  cf  a  lines  sim¬ 
ilar  to  that  suggested  by  Prardt 1  In  block  compression  ((8), 
p.  228)  or  by  Alexander  (l)  in  rolling.  RE  is  assumed  to  be  a 
line  of  stress  discontinuity,  so  that  the  angle  between  the 
sliplines  is  discontinuous  on  it,  while  the  velocity  is  con¬ 
tinuous.  The  suggested  slip  line  is  examined  in  detail  in  Sec¬ 
tion  6. 

E.  Boundary  Conditions 

The  boundary  conditions  or  the  different  portions  of  the 
plastic  flow  domain  ABDETPCA  (Figure  10b)  have  been  in  part  dis 
cussed  in  Section  4,  but  will  be  briefly  recalled  here.  The 
problematic  boundary  condition  -  that  of  stress  on  the  wheel  - 
will  be  discussed  in  more  detail  later. 
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The  separation  line  AB  is  assumed  to  be  an  a  characteristic 
line  with  a  velocity  jump  along  it.  The  line  DE  is  assumed  to 
be  a  horizontal  envelope  of  a  lines  with  a  varying  velocity 
jump  along  it.  The  line  ET  is  an  a  line  with  continuous  veloc¬ 
ity  across,  i.e.,  u  =  1,  v  =  0  on  this  boundary.  CB  and  CD  are 
assumed  to  be  /3  lines  with  continuous  velocity  across  u  =  u^  , 
v  =  0  (Figure  10a).  The  free-surface  RT  is  stress-free  and  a 
streamline  (Equations  [4.l]  and  [4.2],  The  location  of  all  these 
lines  is  a  priori  unknown. 

The  only  remaining  boundary  is  the  soil-wheel  interface 
which  is  obviously  a  streamline.  An  additional  stress  condition 
is  necessary  in  order  to  solve  the  problem.  We  have  not  been 
able  to  find  experimental  evidence  on  this  condition  in  the 
literature  for  the  case  of  a  wheel-clay  interface. 

A  perfectly  rough  wheel  may  be  imagined  as  a  rim  with  a 
series  of  spuds  which  cut  : no  soil.  The  spuds  have  to  be  very 
small  if  the  wheel  is  to  be  considered  as  a  smooth  continuous 
line  and  the  quantity  of  entrained  material  has  to  be  negligible. 
A  perfectly  smooth  wheel  may  be  imagined  as  a  smooth  lubricated 
surface  with  zero  shear  stress  along  it.  A  real  rough  wheel 
is  probably  somewhere  between  the  two.  In  Uffelmann’s  experi¬ 
ments  for  Instance,  (l8),  the  highest  shear  stress  attained  on 
the  wheel  was  t  =  0.3*  In  (5)>  rough  and  smooth  wheels  tested 
in  the  same  kinematical  conditions  gave  approximately  the  same 
resultant  forces.  We  are  inclined,  therefore,  to  admit  that 
the  shear  stress  on  rigid  wheels  is  less  than  the  yield  stress 
and  there  are  no  perfectly  rough  wheels.  In  reality  the  problem 
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is  more  complex.  Between  the  wheel  and  trie  soil  body  moving 
with  different  speeds  a  boundary  layer  of  wet  clay  probably  de¬ 
velops.  This  prob Lem  needs  special  consideration  and  analysis. 

At;  this  stage.,  and  in  the  framework  of  the  rigid-plastic 

model,  it  will  be  assumed  that  on  the  forward  part  of  the  wheel 

AC  tne  srear  stress  T  Is  constant  and  given  (0  *  T  £  l).  Its 

w 

magnitude  has  to  be  determined  experimentally.  It  will  be  shown 
that  there  is  no  major  difficulty  In  solving  the  problem  when 
more  complex  relationships  are  used,  like  ?  depending  on  the 
slippage  along  the  wneel  or  on  pressure;  the  effect  of  such 
relationships  on  tne  plastic  flow  is  left  for  future  investiga¬ 
tion.  At  the  rear  portion  CR  It  is  assumed  tnat  the  shear  stres 

drops  along  the  wheel  from  T  at  C  to  f  =  0  at  R.  This  type  of 

w 

T  distribution  se<-ms  'o  be  confirmed  by  experiments  with  driven 
wheels  (l8) 

In  addition  to  velocity  and  stress  conditions  the  line  ACR 
is  the  only  one  whose  location  Is  known,  having  the  equation 
x3  +  ya  =■  1  The  determination  of  the  silp'llne  field  in  tne  xy 
physical  plane  smarts  always  from  this  line  Tne  angle  of  de- 
tachmert  6^  is  unknown  a^a  constitutes  one  of  ire  Important 
parameters  of  tne  problem 
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6.  THE  SOLUTION  OF  THE  QUASI -STATIC 
EQUATION  FOR  THE  DRIVEN  WHEEL 


A .  The  Method  of  Solution 

The  wheel-soil  problem  belongs  to  the  difficult  category 
of  problems  in  which  both  klnematical  and  stress  boundary  con¬ 
ditions  are  Imposed  on  the  different  boundaries,  which  is  dif¬ 
ferent  from  the  foundation  engineering  problems  in  which  the 
problems  are  formulated  in  terms  of  forces  solely.  For  this 

reason  the  four  variables  p,  0,  v  ,  v0  are  Interrelated  and 

a  p 

the  whole  set  of  Equations  [3.131  and  [3-14]  has  to  be  used. 

The  method  used  here  is  the  analytical  method  of  Hill's 
book  ((8),  Chap.  VI ) .  The  problem  will  be  solved  in  an  auxil¬ 
iary  characteristic  plane  rather  than  in  the  physical  xy  plane. 


Selecting  arbitrarily,  but  conveniently,  two  parameters 

a  and  /3  such  that  a  is  constant  on  /3  lines  and  /3  is  constant  on 

a  lines,  the  variables  p,  0,  v  ,  v. ,  x,y  become  functions  of 

t  ®  P 

a,/3  and  satisfy  the  following  set  of  equations 


|£.  +  2  M  =  0 

9a  S/3 

[6.1] 

^2-  _  2  —  =  0 

9/3  9/3 

[6.2] 

Sva  _  90  n 

9a  V/3  9a 

[6.31 

I* 
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bp  a  bp 


0 


[6.4] 


*2. 

ba 


tan  0  -  0 

aa 


[6.5] 


bp 


+ 


cotan 


0 


[6.6] 


Equations  [6.1]  -  [6.4]  are  obtained  from  Equations  [3. 1,3]  -  C  3  .  l4  3 

by  reducing  by  the  factors  h  =  dx  /do  ,  h.  =  dxD/d£.  h  and  h„ 

ot  cf.  p  p  cl  p 

are  the  metrics  coefficients  (or  Lamme  f  coefficient  s)  and  tney 
represent  the  scale  factors  Detween  the  characteristic  and  phys¬ 
ical  planes.  A  well-known  property  of  Hyperbolic  systems  is 
their  invariance  under  a  cnange  In  characteristic  variables.  Note 
that  tin  full  equations  of  flow  [3.6]  and  [3.7]  do  not  possess 
this  property . 

Equations  [6.5]  and  [6.6]  are  simply  a  transcription  of  the 
characteristic  directions  (Equations  [3. 15]  ) 

A  convenient  transformation  of  Equations  [6,5]  and  [6.6]  is 
achieved  by  using  the  new  variables  ((8),  Chap  VI). 


/  =  x  cos  0  +  y  sin  0  [6.7] 

y  =  -x  sin  0  +  y  cos  0 


which  give 
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*1 

3a 


+  x 


30 

3a 


0 


[6.8] 


[6.9] 


B.  An  Incorrect  Sllpllne  Field 

Before  proceeding  to  the  detailed  solution  of  the  flow 
problem  it  is  worthwhile  to  discuss  a  slipline  field  suggested 
in  the  early  stage  of  this  work  which  proved  later  to  be  incor¬ 
rect.  This  analysis  is  also  valuable  since  in  (22)  it  is  claimed 
that  the  experiments  show  such  a  field. 

In  the  preliminary  stage  of  the  v ork  it  was  assumed  that 
the  plastic  region  is  delimited  by  a  closed  a  line  (Figure  11). 
This  was  proved  ultimately  wrong  for  at  least  two  reasons: 

(i)  At  R  the  free-surface  is  intersected  by  the  CR  a  line 
at  45  degrees  (Figure  11).  According  to  Section  4A  the  region 
TCRT  is  in  a  state  of  compression.  Under  these  circumstances 
the  detailed  solution  of  flow  in  the  triangle  TCRT  (the  method 
is  given  in  Section  6C)  has  shown  that  the  free-surface  is  con¬ 
cave  (Figure  11)  and  the  speed  is  increasing  from  T  to  R  and 
negative  work  is  done  in  the  entire  TCRT  region,  which  contra¬ 
dicts  the  basic  condition  of  plastic  work. 

♦ 

(il)  The  velocity  jump  along  ABCR  is  constant  and  a  typical 
streamline  has  the  shape  depicted  in  Figure  11.  It  is  easy  to 
ascertain  that  along  the  rear  portion  BCR  the  inertial  shear 
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stress  Jump  along  CR  Is  such  that  the  snear  stress  In  the  rigid 
zone  Is  larger  than  the  shear  stress  in  the  plastic  zone,  i.e., 
larger  than  the  yield  stress. 

For  these  two  reasons  this  initial  slipline  field,  although 
esthetically  appealing  at  first  glance,  has  been  abandoned  in 
the  favor  of  that  of  Figure  10b. 

C .  The  Free-Surface  Region  RETR  (Figures  9b  and  11 ) 

Let  us  consider  in  detail  the  free-surface  curvilinear  tri¬ 
angle  TFRT  in  the  physical  plane  (Figure  12a).  A  point  P  has 
the  coordinates  x,y  and  the  characteristic  coordinates  ,  £p . 

The  characteristic  coordinates  are  selected  as  follows:  the 
value  of  £  on  the  a  characteristic  line  P'P  is  0  =  0^,,  l.e.  the 
value  of  0  at  the  Intersection  of  the  a  line  and  the  free-surface; 
the  value  of  a  on  the  £  characteristic  P"P  is  equal  to  0p„  ,  i.e. 
the  value  of  0  at  the  intersection  of  the  £  line  and  the  free- 
surface.  As  long  as  the  free-surface  RT  is  curved  with  no  in¬ 
flexion  points  a  one-to-one  correspondence  exists  between  x,y 
and  a,/S . 

At  point  T  the  free-surface  has  to  be  tangent  to  the  hori¬ 
zontal  unperturbed  rigid  free-surface.  Otherwise  a  velocity 
Jump  6V  will  exist  along  FT  (Figure  11).  Such  a  velocity  Jump 
Induces  an  Inertial  shear-stress  Jump  along  FT  (see  Section  ^C) 
which  Increases  the  shear  stress  in  the  rigid  zone  beyond  the 
yield-stress.  Hence,  the  free-surface  has  to  match  smoothly  the 
horizontal  and  =  0T  =  tt/4. 
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At  poirt  R  the  free-surface  must  be  tangent  to  the  wheel. 

Again,  an  angle  different  of  zero  between  the  free-surface  and 

the  wheel  means  a  discontinuity  in  velocity  along  RF.  But  a 

velocity  Jump  cannot  exist  in  the  interior  of  a  plastic  zone 

since  due  to  inertial  shear  stress  Jump  the  shear  stress  on 

one  side  of  the  line  will  be  different  from  the  yield-stress, 

which  is  absurd.  Consequently  a_  =  /3_  =  0_  =  ir/4  +  0  . 

K  K  rv  r 

The  representation  of  RFTR  in  the  characteristic  plane  a/3 
is  very  simple  (Figure  12b).  By  the  definition  of  a  and  /3  the 
free-surface  is  represented  by  the  line  a  =  /3,  while  the  char¬ 
acteristics  are  parallel  to  the  axes.  Now  p,  0,  v  ,  vfi,  x,  y 

a  p 

as  functions  of  a  and  /3,  for  a  given  0^  have  to  be  determined. 

(i)  Solution  of  p,  0  (dynamics  of  flow)  - 

The  elimination  of  p  from  Equations  [6.1]  and  [6.2] 
shows  that  0  satisfies  the  equation 


d3  0 
3ad/3 


0 


The  general  solution  of  Equation  [6.10]  is 


[6.10] 


0  =  f  (a )  +  h(/3) 


[6.11] 


where  f  and  h  are  arbitrary  functions.  On  the  free-surface 
TR(a  =  /3)0  =  a=  /3by  definition,  which  gives 
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f(/5)  +  h()S)  =  0 
or  h(/3)  =  /3  -  f(/3) 

0  =  f (a)  -  f(jS)  +  jB  .  [6.12] 

From  Equations  [6.1]  and  [6.2]  it  Is  easily  found  that 


P  -  -2 


I 30  .  30 

J  3a  3jS 


d/3 


[6.13] 


The  Equations  [6.12]  and  [6.13]  yield 


p  =  -  2[f(a)  +  f(jB)  -  0]  +  C 


[6.14] 


On  the  free-surface  RT  the  Isotropic  pressure  is  constant  and 
for  the  compression  state  (Equation  [4.3]  p  =  1.  Hence,  for 
a  =  /3  Equation  [6.14]  gives 


-4f(a)  +  2a  +  C  =  1  . 


[6.15] 


From  Equations  [6.12],  [6,l4]  and  [6.15]  the  final  expres¬ 
sions  for  p  and  0  are 
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0  = 


a  +  ft 
2 


p  =  -a  +  ft  +  1 


[6.16] 

[6.173 


which  completely  solves  the  problem. 

It  is  worthwhile  to  note  that  at  point  F  (Figure  12) 

ap  +  n 

0p  =  - g ^  +  —  >  i.e.,  the  characteristic  line  TF  bends 

0 

with  an  angle  —  and  the  characteristic  FR  with  the  same  angle 
0 

—  .  The  pressure  drops  along  TF  and  reaches  a  minimum  value 

at  F,  p„  =  -a  +  ft_  +  1  =  1-  0  which  means  that  the  soil  is 
r  r  r  r 

under  slight  tension  in  this  region. 

(li)  Solution  of  v  ,  v.  (kinematics)  - 

u  ft 

From  Equations  L 6 .3 3 ,  [6.4]  and  [6.l6]  it  is  found  that 
v^  and  v^  satisfy  the  following  relationships 


dv  v0 

-SL  _  J.  =  o 

3a  2 


L  6 . 18] 


3v.  v 

— ^  +  —  =  0 
aft  +  2 


or,  by  elimination  of  v  and  vQ 

a  p 
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aav  v 
a  ,  a 

+  T  =  0 

&a  v  v 

_Jk  +  o 

-v  -s  ^  '  1.  ^ 


[6.191 


l.e.,  v  and  v.  satisfy  the  telegraph  equation. 

cip 

The  kinematic  boundary  conditions  are  as  following:  On  the 
characteristic  FT  (Figure  12)  the  velocity  is  continuous  and 
equal  to  the  rigid  velocity  u  =  1.  Hence,  from  Equation  [3.12] 


7T 


+  a 


v  =  cos  0  =  cos 
a 


sin  0  =  -  sin 


7T 


+  a 


[6.20] 


On  the  free-surface,  according  to  Equation  [4.4] 


v  +  vfl  =  0  .  [6.21] 

d  p 

v  and  v0  may  be  found  from  Equations  [6.19],  16.20]  and 
ex  p 

[6.21]  by  the  following  procedure.  Let  us  introduce  a  new  var¬ 
iable  v  =  v  +  v„  .  v  satisfies  the  following  equations  and 
ex  p 

boundary  conditions  in  the  characteristic  plane  (derived  from 
Equations  [6.19],  [6.20]  and  [6.21] 
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aav 

a  £ 


+ 


V  =  cos 


v  =  0 


^I< 

II 

0 

(in  TFRT) 

[6.22] 

TT  ,  7 r 

IT  +  a  r  +  a 

2  Sin  2 

(on  TF) 

[6.23] 

.  (on  TR) 

[6.24] 

The  triangular  domain  TFRT  in  the  a,£  plane  (Figure  12b)  may  be 
reflected  across  TR  in  TF'RT.  From  the  symmetry  of  the  tele¬ 
graph  equation  and  the  condition  [6.24]  v  may  be  determined  in 
the  whole  rectangular  domain  TFRF'T  by  reflection.  Since  TF'  is 
the  reflection  of  FR,  v  takes  on  the  value 


v 


cos 

i 


sin 


(on  TF ' ) 


[6.25] 


The  solution  of  the  telegraphic  equation  for  a  function  given 
along  two  intersecting  characteristics  by  the  Riemann  procedure 
is  described  by  Hill  ((8),  p.  155).  In  the  present  case  the 
solution  is  (Figure  12b) 


v(P) 


[6.26] 


where  G(a,£ ; 


a 

P 


~\J (ap  -  °0(£p  -  0) 


is  the  Green  Function. 
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The  substitution  of  the  expressions  of  the  Green  function 
and  the  values  of  3v/£a  and  dv/d|S  on  the  boundaries  in  Equa¬ 
tion  [6.26]  gives  finally 


Equation  [6.2?]  completely  determines  v.  In  order  to  find 
for  Instance,  the  second  equation  of  [6.18]  is  written  as 


3j3  2  2 


[6.28] 


The  integration  of  this  ordinary  differential  equation  along  a 
j3  line  from  ft  to  P  (Figure  12)  gives 

J  +  a  2  ^  “  ¥'  j3/2  f  -5/2 

v^(a,j3)  =  -  sin  — ^ —  e  -  c  J  v(a,S)  e  d?  , 

tt/4  L  6 . 29l 

which  determines  v.  once  v  is  known.  Obviously  v  is  found  from 

p  a 
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The  integrals  in  Equations  [6.271  and  [6.29]  have  been  com¬ 
puted  numerically  (by  Simpson  rule)  on  HYDRONAUTICS •  IBM  1130 
computer . 

The  resulting  values  of  v  ,  v_,  0  and  V  =  V v  3  +  v_a  are 

a  p  a  p 

given  in  Appendix  1.  The  computations  have  been  carried  out  to 

the  maximum  possible  value  of  8^  =  tt/4  and  the  results  may  be 

used  for  any  arbitrary  0^  in  the  range  0  <  0^  <  tt/4.  Each  group 

of  results  in  Appendix  1  corresponds  to  a  constant  a  in  the 

range  ir/ 4  <  a  <  ir/2  and  a  variable  J3  in  the  range  tt/4  *  /3  <  a. 

The  velocity  V  on  the  free-surface  has  the  value  of  the  last 

line  of  each  group,  where  a  =  .  It  is  easy  to  ascertain  that 

the  velocity  decreases  along  the  free-surface,  or  increases  from 

T  (V  =  l)  to  R  (V_  =  1.75  for  0  =  tt/4).  For  moderate  values 

1  k  r 

of  0^,  as  encountered  in  applications,  VR  is  smaller.  The  ve¬ 
locity  distribution  along  TR  will  be  used  in  Section  7A  in 
order  to  determine  the  magnitude  of  the  recovery  angle  0^. 

The  results  of  this  section  and  of  Appendix  1  are  valid  not 
only  in  the  triangle  RFTR,  but  in  the  whole  region  REFTR  (Fig¬ 
ure  10b).  The  plastic  flow  in  RFER  is  just  the  continuation  of 
the  flow  in  RFTR  across  the  characteristic  RF. 

(iii)  The  determination  of  x,y  (geometry  of  flow)  - 

The  mapping  of  the  free-surface  triangle  TFRT  from  the 
characteristic  plane  to  the  physical  plane  may  be  done  by  using 
Equations  [6.5],  [6.6]  or  Equations  [6.8],  [6.9].  In  order  to 
solve  the  problems  the  location  of  the  line  RE  (Figure  10b)  has 
to  be  known.  In  principle  this  can  be  done  by  starting  with  the 
solution  of  the  x,y  equations  from  the  wheel,  but  this  is  not 
done  here. 
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The  x(a,0)  and  y(a,/3 )  of  the  line  RE  are  enough  in  order  to 
solve  the  problem  since  an  additional  boundary  condition  exists 
on  the  free-surface,  namely  the  fact  that  0  -  0  =  tt/4  (Equa¬ 
tion  [4.31.  Since  tan  0  =  dy/dx,  on  the  free-surface  RT 


=  tan 
dx 


tan 


a  +  ft 
2 


[6.30] 


This  additional  condition  permits,  in  principle,  the  mapping 
of  the  whole  region  TERT  from  the  a,/3  plane  onto  the  x,y  plane. 

D.  The  Bow  Region  ABCA  (Figures  10b  and  13) 

In  this  case  the  characteristic  coordinates  of  a  generic 
point  P  (Figure  13a)  are  selected  as  following:  a  is  the  angle 
„  at  the  intersection  between  the  /3  line  PP"  and  the  a  base 
line  AP"B  and  /3  is  0p(  at  the  intersection  of  the  a  line  PP'  and 
the  /3  base  line  BP'C. 

Some  general  relationships  may  be  found  before  solving  the 
problem  in  detail.  At  point  A  (the  leading  edge)  the  velocity 
inside  the  plastic  zone  has  to  be  parallel  to  the  wheel,  since 
no  velocity  discontinuity  is  admitted  inside  the  plastic  zone 
from  inertial  considerations  (see  Section  4d).  AB  is  a  line  of 
velocity  Jump.  The  Jump  5V  may  be  found  from  continuity  con¬ 
siderations  at  point  A  (Figure  13b). 

From  the  velocity  triangle,  by  projecting  on  the  direction 
of  6V  and  on  the  normal  to  it,  the  following  is  obtained: 
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sin  0 
~  sin  (0^ 

sin  0 

BV  =  Vo  COS  (0.  -  0  )  -  COS  0.  =  -r-rx - — — Q — r  [6.32] 

a  '  l  a'  I  sin(0.  -  0  ) 

8l 

The  angle  0  is  related  to  the  sinkage  z  by  (Figure  12a) 

8l 

Z  =  1  -  COS  0  .  [6.33] 

fit 

The  angle  0^  is  determined  by  the  shearing  stress  on  the 
soil-wheel  Interface  at  A .  0^  is  the  angle  between  the  a  line 

AB  and  horizontal,  i.e.,  0^  =  0 wA<  The  relationship  between 
and  0wA  is,  according  to  Equation  [4.71 

=  cos  2(0w  -  9)  =  cos  2(0^  -  0  )  [6.34] 

which  determines  0^  for  given  and  0&.  Obviously  may  vary 

in  the  range  0  *  Tw  s  i.e.,  between  perfectly  smooth  and 

perfectly  rough  conditions.  For  =  0  at  A,  0^  -  0&  =  tt/4 

while  for  T  =1,  0,  -0  =0  and  the  characteristic  AB  sticks 

w  l  a 

to  the  wheel. 

According  to  our  present  concept  on  soil-wheel  interface 

(see  Section  5E)  T  is  assumed  constant  on  AC. 

w 


[6.31] 
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The  line  BC  is  a  characteristic  £  line  between  the  plastic 
zone  ABCA  and  the  assumed  rigid  core  BCDB  (Figure  10b).  The 
velocity  in  the  rigid  core  is  horizontal  and  constant  u  =  u  , 
v  =  0.  There  is  no  velocity  jump  along  BC  since  at  C,  inside 
the  plastic  zone,  the  velocity  along  the  wheel  is  also  hori¬ 
zontal.  Hence,  the  slipline  AB  has  a  horizontal  tangent  at  B 
(0_.  =  0).  Since  &V,  the  velocity  jump  along  AB  is  constant 

Jd 

(Section  4c),  we  have 


ut>  °  1  +  BV  *  1  +  alrffi/  -  8  )  ’  t6-351 

v  Si 

which  gives  ufe  in  terms  of  known  quantities.  From  continuity 
considerations 


or 


z  +  h  =  u^h 

z  (1  -  cos  9a)  sln  *  ®a ) 

h  =  u  -  1  =  sirTe 

b  a 


[6.38] 

C  6.373 


which  gives  the  thickness  of  the  plastic  zone  at  bottom. 

Finally,  for  the  slipping  wheel  considered  here  the  wheel 

peripheral  velocity  V  1  =  uu'r'  has  to  be  larger  than  the  soil  ve- 

w  . 

locity  at  the  interface  AC  (Section  5E)  •  Since  the  maximum  soil 
velocity  occurs  at  bottom  where  u  =  u^,  the  following  inequality 
holds 


V  £ 
w 


V  - 
w 


V 

w 


1 


or  in  terms  of  slippage  coefficient  S  = 
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Using  the  value  of  of  Equation  [6.351 

sin  0 

£L 

^  ^  sin(0 .  -  0  )  +  sin  6 

v  Q-  a. 


[6.381 


[6.391 


This  is  the  minimum  slippage  necessary  in  order  to  ensure 

a  positive  T  along  the  wheel.  The  real  slippage  may  be  larger 
w 

than  and  rigid-plastic  considerations  solely  permit  to 
remain  unchanged.  A  more  detailed  study,  experimental  or  theo¬ 
retical,  may  reveal  a  relationship  between  and  T^. 

If  we  assume  that  the  real  slippage  coefficient  S  is  equal 

to  S  ,  then  a  relationship  between  T  and  S  (for  given  z)  is 
m  r  w  m 

established  through  Equation  [6.391-  For  T  “  1  >  S  =  1,  as 

w  m 

expected . 

Now,  a  detailed  solution  of  the  plastic  flow  is  presented 
along  the  lines  of  Section  6D.  The  domain  ABCA  is  represented 
in  the  characteristic  a,£  plane  in  Figure  13C.  The  representa¬ 
tion  follows  straightforward  from  the  definition  of  a  and  j3 .  The 
wheel-soil  interface  AC  is  represented  by  an  unknown  curvilinear 
line  in  the  characteristic  plane. 
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(i)  The  determination  of  p,  0  (dynamics) 

Again,  0  has  the  general  expression  C6.ll],  But,  from 
the  definition  of  a  and  £,  we  have  0  =  a  on  AB  (0  =  0)  and  0  - 
on  BC  (a  =  0).  This  gives  immediately 

0  =  a  +  j9  [6.40] 

The  isotropic  pressure  p  is  related  to  0  through  Equation 
[6.13].  Hence, 


p  =  -2(a  -  £ )  +  const  [6.4l] 

Let  pc  be  the  pressure  at  C.  Then,  [6.4l]  becomes 

p  =  p  -  2(a  -  fi  -  0  )  [6.42] 

C  wc 

where  0  =0  -8=_A  cos-1  T  . 

wc  w  *  w 

The  pressure  increases  along  the  wheel  from  p~  to  p  ,  p 

v/  A  A 

being  given  by 

PA  =  ?o  +  2(6t  +  '  C6-'131 

The  normal  stress  acting  on  the  wheel  is  (Equation  [3.7]) 

o  =  -  p  +  sin  2(0  -6). 

w  w 


[6.44] 
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Slnce  T  =  cog  2(0  -  9 )  is  assumed  constant,  0  -  0  is  constant 

w  x  w  '  w 

and  ow  is  simply  equal  to  the  isotropic  pressure  plus  a  constant. 

For  a  perfectly  smooth  wheel  (T  =  0),  0  -  0  =  -ir/4  and 

'  w  '  w 

o  =  -p  -  1;  for  a  ~ough  wheel  (?  =  1 ) ,  0  -0=0  and  o  =  -p . 

w  v  w  w  w 

(ii)  The  determination  of  v^,  (kinematics)  - 

Similarly  to  Equations  [6.18]  and  [6.10]  we  have  this 
time  for  v  and  v„ 


dv 

“  v  0 

da  p  “ 

dvfl 

°  +  V  =0 
dp  a 

C6.U5] 

d3v 

a 

STS'  +  V  =0 
dadp  a 

*3VP 

dadp  +  Vp  ~  °  * 

[6.1160 

The  boundary  conditions  for  the  telegraph  Equations  [6.461 
are  as  follows: 

-  On  AB  and  satisfy  the  conditions  along  a  line  of 
separation  between  a  plastic  and  a  rigid  zone  (Section  4C) 


v  =  1  +  5V  cos  0 

“  t.6.117] 

v .  =  -  sin  0 
P 


or  by  Equation  [6.4o] 
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v  =  1  +  6V  cos  a 
a 

(on  /3  =  0;  -0^  <  a  <  0)  [6.483 

v  =  -  sin  a  . 

P 


-  On  the  line  BC  there  Is  no  velocity  Jump,  hence 


v  =  u,  cos  0  »  u,  cos  6 
a  b  b 


-u,  sir.  0  =  -u.  sin  8  . 

D  D 


(on  a  =  0;  -0  ^  <  p  <  0) 
w  c 


[6.493 


These  conditions  completely  determine  the  values  of  vq  and 

v  at  any  point  in  the  characteristic  rectangle  ABCB'A  of  Flg- 
P 

ure  13c). 


Following  the  same  Riemann  method  as  In  Section  6d,  v^  and 

vfl  at  a  point  P(a,p)  are  given  by 
P 


P" 


B 


da 


B 


dv 

G  jf  dj3  r  G(P,B)  va(B) 


P' 


yp) 


dv 

— ■  dp  +  G(P,B)  v^(B) 


P' 


L  6 . 30  J 
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where  G,  the  Green  Function,  is 


0(ap,|9p;o,/J)  -  J0  2  V(<ya)(/>p-/S) 


E  6.51] 


Using  the  values  of  v  ,  v.  on  the  boundaries  (Equations  [6.483 

ct  p 

and  [6.49]  and  the  definition  of  G  (Equation  [6.5l3)»  Equations 
[6.503  become 


V 


a 

r  r 

0 

B  * 

•  /  J0 

sin  SdS  -  ub  /  jl 

2  Va(fS-S) 

sin  5d? 

J  l 

0 

^  JQ  (2^) 

?  J  * 

0 

[6.523 

ra  r 

-  J  2V(a-S)/S 

cos  5d5  -  ub  I  JD 

2Va(^-5) 

cos  5d? 

J  1  J 

0 

J 

0 

[6.533 

For  the  illustration  of  the  method  we  have  selected  an  ex¬ 
ample  with  the  following  data 


9  «  30  -  z  =  0.134 

a 


0-0=0  =  -35°  -  T  =  0.34  ,  0.  *  65°  [6.543 

w  wc  w  7  l 


v  and  V-  have  been  determined  by  numerical  integration  and  their 

ot  P 

values  for  different  a  and  /3  in  the  rectangle  -0^  <  a  <  0  , 

0  <  £  <  0  are  given  in  Appendix  2. 

w  c 
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The  location  of  the  line  AC  is  found  now  by  using  the  con¬ 
dition  that  AC  is  a  streamline  and  a  line  of  constant  T  *  T  . 

w 

According  to  Equation  [4.8],  v  and  v0  satisfy  the  condition 

a  p 

vA 

-  tan  (0  -  0  )  [6.553 

va 

or,  in  our  example  (0  -  0  =  -35°) 

w 

v 

^  -  tan  35°  =  0.7  .  [6.56] 

a 

In  the  numerical  results  v./v  has  been  given  for  each  point 

p  a 

(Appendix  2)  and  it  was  a  simple  matter  to  pick  from  the  results 
the  location  of  the  points  which  satisfy  Equation  [6.56].  In 
Figure  l4  the  characteristic  plane,  including  AC,  is  represented 
according  to  the  data  of  Appendix  2. 

Once  AC  is  known  the  same  Appendix  2  provides  the  speed 

v  ="^7*-  +  vQ3  along  the  wheeL.  In  Figure  15  the  velocity  varla- 
cx  p 

tlon  along  the  wheel  and  along  corresponding  points  (l.e.,  on 
same  a  lines)  along  the  separation  line  AB  is  given  graphically. 
In  Figure  17  the  velocity  distribution  along  a  few  character¬ 
istics  is  also  given. 

The  velocity  along  the  wheel  increases  from  its  value  at 
the  leading  edge  A  (V  =  1.57)  to  its  value  at  the  bottom  B 
(V  =  1.88)  and  the  requirement  of  positive  work  is  satisfied. 
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The  velocity  along  the  slipline  AB  is  somehow  higher  than  the 
value  on  AC.  A  rough  check  has  shown  that  the  positive  work  con 
dition  (Equation  13.22])  is  satisfied  everywhere  in  the  plastic 
zone.  The  small  difference  between  V  on  AC  and  AB  for  the  same 
£  slipline  shows  that  the  distribution  is  almost  uniform  on  the 
j3  lines,  but  changes  in  a  direction  (Figure  17). 

(ill)  The  solution  of  x,y  (geometry)  - 

Equations  [6.8]  and  [6.9]  have  been  used  in  order  to 
map  the  characteristic  region  of  Figure  15  onto  the  physical 
plane. 

Along  the  wheel  AC  x  and  y  are  known,  namely, 
x  =  sin  0  y  =  -  cos  0 

According  to  Equations  [6.7]  x  and  y  are  given  by 

x  =  sin  0  cos  0  -  cos  0  sin  0  =  -  sin  (0  -  0 ) 
w  w  w 

y  =  -  sin  0  sin  0  -  cos  0  cos  0  =  -  cos  (0  -  0 ), 

w  w  w  ' 

But,  from  the  assumption  of  constant  we  have 

0-0=0  =  -  \  cos"1  T  [6.591 

w  wc  2  w 

Hence,  from  Equation  [6.5*0 

x  =  sin  35°  y  =  -cos  35° 


[6.573 


[6.58] 


% 


(On  AC) 
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With  these  values  of  x  and  y  along  AC,  the  Equations  [6.8] 
and  [6.9]  have  been  integrated  numerically,  by  finite  differ¬ 
ences.  The  results,  in  terms  of  the  mapping  of  a,£  region  ABCA 
on  x,y  plane,  are  given  in  Appendix  3.  The  slipline  field  has 
been  plotted  graphically  in  Figure  l6.  The  sliplines  are  ortho¬ 
gonal,  as  required  by  the  theory.  They  intersect  the  wheel  at 
the  constant  angle  corresponding  to  =  0.34  ( 0 -  0  =  -35°). 

Hence,  for  the  selected  example,  the  plastic  flow  in  the 
bow  region  has  been  completely  solved.  The  stress  distribution 
has  been  determined  (Equations  [6.40]  and  [6.42]  excepting  the 
constant  p  which  has  to  be  found  from  the  matching  of  bow  and 

c 

rear  plastic  regions. 

There  is  no  problem,  in  principle,  to  find  graphically  the 
streamlines  in  the  bow  region,  but  the  work  involved  is  quite 
tedious.  A  simpler  problem  is  to  depict  the  velocity  distribu¬ 
tion  along  characteristics,  based  on  the  data  of  Appendices  2 
and  3.  Such  a  representation,  for  the  discussed  example,  is 

given  in  Figure  17.  Figure  17  shows  that  v  is  almost  uniformly 

c& 

distributed  along  the  j3  characteristics,  while  vD  has  approxi- 

P 

mately  a  linear  variation.  This  observation  may  help  in  con¬ 
structing  the  kinematical  field  approximately  without  solving 
the  problem  in  detail.  It  may  be  also  useful  in  future  work  in 
which  inertial  terms  inside  the  plastic  zone  or  work-hardening 
effects  would  be  taken  into  account. 
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E.  The  Plastic  Zone  CDERC  and  the  Rigid  Core  CBDC  (Figures  10b 
and  18) 

The  solution  of  the  plastic  flow  in  CDERC  has  to  start  with 
the  horizontal  envelope  of  a  lines  DE.  There  is  no  systematical 
study  of  such  lines  in  literature. 

Again,  the  representation  of  the  plastic  domain  in  the  char¬ 
acteristic  plane  has  to  start  with  a  selection  of  the  character¬ 
istic  variables  a,jS.  The  zone  CDERC  is  represented  in  detail  in 
Figure  l8a. 

Let  a,jS  in  the  subregion  DQED  be  the  following  (Figure  l8a): 

a  is  the  distance  between  D  and  P"  on  the  base  line  DF  (a  =  s  „ 
P  P  P 

and  jSp  is  Sp,.  A  one-to-one  correspondence  is  thus  established 

between  points  in  DQED  and  the  a,£  plane.  The  image  of  DE  in 

a,jS  plane  is  obviously  the  line  a  =  jS  (Figure  l8b). 

Since  0  =  0  on  this  line  it  is  easy  to  ascertain  that  0  has 
the  expression 


0  =  a  -  jS 


L  6 . 6o] 


in  the  domain  DQED. 

The  Equations  [6.3],  [6.4],  [6.8]  and  [6.9]  become  now 
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0 


a/3 


v  =  0 
a 


C6.6l] 


3a 


+  x 


0 


3x 

a/3 


+ 


y 


o 


All  the  functions  v  ,  v_,  x  and  y  satisfy  telegraphic  equa- 

a  p 

tlons  of  the  type 


33  f 
3ad/3 


f  =  0 


[6.62] 


which  admits  as  Green  function 


G(a,0;ap  f/3p  )  =  Iq  I  2  V(ap-a}(^  -^  ) 


[6.63] 


Hence,  all  the  dependent  functions  may  be  found,  in  principle, 
by  Riemann's  method. 

The  plastic  domain  CDQR  (Figure  l8a)  may  be  represented  on 
the  same  a,£  plane  by  selecting  as  a  variable  the  0  angle  on  the 
base  line  CD.  With  0  =  a  -  0  the  region  CDQRC  is  mapped  on  the 
lower  part  of  the  a,/3  plane  (Figure  l8b ) . 


HYDRONAUTICS,  Incorporated 


-69- 

Both  images  of  CR  (the  wheel)  and  o.f  RQE  (the  stress  dis¬ 
continuity  line)  in  the  a,/S  plane  are  unknown,  including  the  de¬ 
tachment  angle  0  . 

The  functions  v  ,  vc ,  x,y  and  the  location  of  CR  and  RE  in 
a  p 

the  characteristic  plane  have  now  to  be  found  by  using  Equa¬ 
tions  [6.6ll  and  the  boundary  conditions.  The  latter  are  the 
following:  (i)  On  CD  the  velocity  is  known  (u  =  ufe,  v  =  0^. 

(li)  On  DE  one  of  the  velocity  components  is  known,  v  =  v„  -  0. 

P 

x  and  y  are  also  known,  x  =  s,  y  =  0.  Hence,  since  0  =  s,  one 
has  x  =  a3  +  £3  =  a  J~2,  y  =  0.  (iii)  On  RE  the  velocity  com¬ 
ponents  v  ,  v.  have  to  match  continuously  the  velocity  in  the 

CL  p 

free-surface  region  RETR  which  have  been  found  previously.  More¬ 
over,  since  RE  is  a  line  of  stress-discontinuity,  it  bisects  the 
angle  between  the  sliplines  on  its  two  sides  (Figure  l6a).  This 
is  an  additional  condition  of  matching  between  the  two  plastic 
zones.  (lv)  The  wheel  CR  is  a  streamline  and  x,y  are  interre¬ 
lated  by  the  Equation  x3  +  y3  =1. 

Unfortunately  the  analytical  solution  of  the  problem  by 
Rlemann's  method  is  very  difficult  since  the  boundary  conditions, 
unlike  the  bow  and  free-surface  regions,  do  not:  permit  to  march 
with  the  solution  from  one  boundary.  The  analytical  approach 
leads  to  integral  equations  for  the  dependent  functions  which 
have  to  be  solved  by  matching  the  solution  with  the  flow  in  the 
free-surface  region. 

A  finite  difference  solution  is  also  difficult  because  of 
the  unKnown  location  of  CR  ana  RE.  For  this  reason  a  detailed 
solution  of  the  plastic  flow  in  the  region  CDERD  has  been  abandoned 
at  this  stage.  An  approximate  solution  is  presented  in  Section  7. 
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Finally,  the  geometry  of  the  rigid  core  CBDC  is  determined 
from  the  solution  of  the  flow  in  the  bow  region  from  one  side 
(CB)  and  from  the  solution  of  the  flow  in  CREDR  from  the  other 
side  (CD). 

The  stresses  on  CR  are  found  from  the  solution  in  CREDC 
starting  from  the  free-surface  where  p  =  1.  The  pressure  on  the 
left  hand  side  of  the  rigid  core  at  C  has  to  be  found  from  the 
equilibrium  of  the  rigid  core.  When  pc  is  known,  the  stress  dis¬ 
tribution  in  the  bow  region  is  totally  determined. 
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7.  THE  APPROXIMATE  DETERMINATION  OF  THE 

DETACHMENT  ANGLE  0  AND  OF  STRESS 

r 

DISTRIBUTION  ALONG  THE  WHEEL 

The  detailed  solution  of  the  problem  of  plastic-flow  beneath 
a  wheel,  discussed  and  exemplified  in  section  6  requires  a  large 
amount  of  computations.  Since  the  solution  is  an  approxi¬ 
mate  one  because  of  the  different  simplifying  assumptions  under¬ 
lying  the  mathematical  representation,  it  is  highly  desirable  to 
find  a  simple,  approximate  method  of  determining  the  forces  on 
the  wheel.  Such  a  method  is  presented  in  this  section. 

A.  The  detachment  (or  recovery)  angle  0  ^ 

6  ^  may  be  found  approximately  from  the  continuity  require¬ 
ments  between  the  0  line  CD  and  the  line  RE  (Figures  10b  and  19). 

Since  there  is  no  soil  flow  across  CR  (the  wheel)  and  DE, 
the  flux  through  CD  and  RE  have  to  be  equal,  i.e. 

R 

h  uw  V  cos  T  ds  [7  1] 

bJ 

E 

Since  the  mapping  of  RE  on  the  physical  plane  is  not  known,  the 
relationship  [7.1]  in  its  exact  form  is  of  little  help  An 
approximate  evaluation  of  the  integral  may  be  obtained  as  follows. 
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First,  Equation  [7*1]  is  equivalent  to 

R  hR 

uvh  =  /V  cosqdsr/  V  dy 


[7.2] 


q  being  the  angle  between  the  normal  to  RE  and  the  velocity 
vector  and  (-  the  angle  between  the  normal  and  x  axis  (Fig.  19). 

At  R  it  is  reasonable  to  assume  that  the  stress  jump  is 

negligible  and  RE  is  tangent  to  the  characteristic  RF.  q  at  R 

is  l..,own  since  the  characteristic  RF  intersects  the  free  surface 

at  tt/4.  Hence  £_  =  7r/4  +  0  ,  while  q„  =  9  .  From  the  other  hand 

V  at  R  is  the  speed  along  the  free  surface  and  it  has  been  found 

by  solving  the  kinematic  problem  in  RTFR  (Section  6c).  For  each 

given  9  ,  V  may  be  found  in  Appendix  1  (in  the  last  line  of 
r  n 

each  group  of  results)  taking  into  account  that  =  tt/4  +  9 at  R. 

H  .>nce  the  integrand  of  [7.2]  has  a  well  determinated  value  at  R, 

for  a  given  0  . 

r 

At  E  again  V  =  u  =  1  is  known  and  there  also  q  =  (•  since 
the  velocity  vector  is  horizontal.  The  value  of  V(cos  q/cos  I* ) 
being  known  at  the  two  extremities,  4  reasonable  approximation 
of  the  integral  [7.2]  is  obtained  by  assuming  a  linear  variation 
of  the  integrand  with  y.  Hence,  Equation  [7.2]  becomes 


1  + 


V_.  cos  9 
R  r 


h 


R 


cos  (tt/4  +  6  ) 


[7.3] 
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A  few  values  of  VD  as 

K 

a  function 

of  0  ,  from  Appendix  1 

are  given  in  the  following 

table 

0 

4> 

V„ 

r 

R 

R 

5.1° 

50.1° 

1.08 

11.2° 

56.25° 

1.20 

20.2° 

65.20° 

1.35 

27° 

72° 

1.76 

31.5° 

76.5° 

1.54 

The  quantity  hR  is  simply  related  to  h  (Pig.  19)  by 


h„  =  h  +  1  -  cos  0 
R  r 


[7.4] 


Equations  [6.35]  and  [6.37]  relate  u,h  to  0  and  0,  -  0  =  -<t> 

b  a  l  a  wc 

Hence,  in  Equation  [7.3]  all  quantities  are  functions  of  0,0, 

r  3 

and  <t>wp  and  it  permits  the  determination  of  0^  (recovery  angle) 
for  given  0g  (sinkage)  and  <t> (shear  stress  at  bow).  Equa¬ 
tion  [7.3]  has  been  solved  by  trig  and  error.  In  Figure  20 
the  dependence  of  0^  on  the  shear  stress  on  the  wheel  is  given 
for  three  sinkages  z  =  O.O36,  z  =  0.134  and  z  =  0.293  corresponding 
to  0  =  15°i  0  =  30°  and  0  =  45°  respectively.  The  minimum 

slippage  as  function  of  z  and  tw  is  given  in  Figure  23. 

The  examination  of  Figure  20  shows  that  0^  increases  with 
sinkage  and  with  t  .  For  small  sinkages,  as  encountered  in  ap¬ 
plication,  and  t  «  0.5,6  is  of  order  0  /2. 

1  w  ’  r  a' 
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B.  The  Stresses  Along  RC  (Figure  10b) 

As  stated  in  Section  the  shear  stress  is  assumed  constant 

r  =  along  the  bow  region  AC,  and  dropping  linearly  along  RC 

from  t  =  t  to  t  =  0 
w 

T  =  Tw  (i  -  e/er)  (on  RC,  o  <  e  <  0r)  [7.5] 

The  isotropic  pressure  at  C,  on  the  right  side  of  CD  (Fig¬ 
ures  12b  and  l8a )  is  found  from  the  integration  of  Equations  [6.1] 
and  [6.2]  along  the  sliplines  RN  and  NC  as 


1  +  2 


*  ) 
wc ' 


[7.6] 


the  terms  in  brackets  being  the  variation  of  <t  between  R-N  and 

N-C  respectively.  The  value  of  is  unknown,  but  its  range  of 

variation  Is  0  <  <t>  .  <  <t>  ,  corresponding  to  N  —  D  or  to  N  -*•  C, 

N  wc 

respectl  /ely.  Assuming  the  average  value  <t>N  =  <J)wc/2  we  obtain 
from  Equation  [7.6] 


P0  =  1  +  tt/2  +  29 
Cr  r 


[7.7] 


The  normal  stress  on  the  wheel  at  C,  on  the  right  side  of 
CD,  is  [Equation  3-7] 
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whero 


Cr 


=  -  p „  +  sin  24 >  =  -(1  +  J  +  20  )+  sin  2<t> 

HCr  wc  2  r  wc 

related  to  t  by  4>  =  -  1/2  cos  1r  . 

w  wc  '  w 


[7.8] 


The  type  of  variation  of  a  from  a  =  0  at  R  to  o  =  a_  at 
17  K  w  w  w  Cr 

C  is  not  known.  From  the  assumed  slipline  distribution  it  is 
reasonable  to  assume  a  linear  variation 


a 


w 


1 


1  +  f  +  20 

2  r 


sin  2<t> 


wc 


(on  CR,  0  <  0  <  0  ) 


[7.9] 


Equations  [7.5]  and  [7.9]  Gave  the  approximate  distribution 

of  shear  and  normal  stresses  along  the  rear  portion  of  the 

wheel  CR  as  a  function  of  6,  0  ,  and  t  (or  equivalently  <t>  ). 

r  w  '  wc 

C .  The  Equilibrium  of  the  Rigid  Core  CBDC  (Figure  10b) 

In  the  rigid  (elastic)  core  the  maximum  shearing  stress  has 
to  be  lesser  than  the  yield  stress  anywhere.  The  shape  of  the 
two  lines  CB  and  CD  is  only  partially  known.  CB  is  known  from 
the  solution  of  the  bow  region  (Figure  l6)  while  CD  has  to  be 
found  from  the  solution  of  the  rear  region.  Figure  16  shows 
that  an  acceptable  approximation  is  to  replace  the  lines  CD  and 
CB  by  two  arcs  of  circles  (Figure  21).  The  two  arcs  are  com¬ 
pletely  determined  by  the  known  angles  at  C(o  )  and  at  B  and  D 

w  0 
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(tt/2)  and  by  the  height  h.  The  isotropic  pressure  distribution 
is  also  determined  from  the  /3  characteristic  relationship  [6.2] 
as  (Figure  19) 


p  =  PCr  ^  20wc  “  2r]  ^on  CB^ 

p  “  PC£  "  2*wc  +  2t>  {on  CD) 


[7.10] 


The  shearing  stress  distribution  along  the  three  sliplines  CB, 
BD  and  DC  is  also  known  and  given  in  Figure  21. 


All  the  stresses  of  Figure  21  are  known  excepting  p^,  the 
pressure  at  C  on  the  left  side  of  CB.  The  equilibrium  of  the 
forces  acting  on  the  rigid  case  in  horizontal  direction  will 


[7-10] 


The  same  type  of  computation  gives  for  the  pressure  resultant 
on  CB 


1  cos  $ 

=  h  p_.  +  2  — — — 
CB  Cl  sin  <t> 


_ 2 _ 

sin  $ 


t7.ll] 
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The  resultant  of  the  shearing  stresses  Is 


cos  4> 


H  =  Hi  =  4h  - - - - 

s  sin  <t> 


wc 


wc 


The  equilibrium  requires 


[7.12] 


HCD  +  Hs  -  HCB  -  0  [7-131 

Introducing  the  expressions  [7.10],  [7.H]  and  [7.12]  Into  [7-13] 
and  after  some  algebraic  manipulations,  it  is  found  that 


PC l  PCr  PC  [7.14] 

i.e.  the  pressure  at  C  is  transmitted  unchanged  though  the  rigid 
core.  This  shows  that  the  rigid  core  is  effectively  not  in  a 
plastic  state,  in  which  case  p  would  increase.  Hence  from  Equa¬ 
tions  [7.7}  and  [7.14] 


P,  -1+1+  20r  17.15] 

D.  The  Stresses  Along  CA  (Figures  10b  and  13a) 

A  method  of  determining  the  stress  distribution  by  the  exact 
integration  of  the  quasi-static  equations  has  been  presented  in 
Section  6d.  The  inspection  of  Figure  14,  which  shows  the  image 
of  the  line  AC  in  the  characteristic  plane,  reveals  that  the 
curve  may  be  approximately  replaced  by  a  straight  segment  between 
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A  and  C.  The  equation  of  this  straightline  Is  (Figure  13c) 


P  +  =  0 


[7.16] 


From  the  other  hand  <t>,  p  and  0  are  related  to  a,fi  by  Equa¬ 
tions  [6.40],  [6.42]  and  [6.58].  From  these  equations  and  [ 7 . 1 6 ] 
we  find  that  along  CA 


a  -  <t> 


2  1  +  —  a 

6  f 

&  1 


[7.17] 


Replacing  p  by  Its  value  of  Equation  [  7 . 1[> ] ,  p  Is  ex- 
0 

pressed  as  function  of  9  along  the  wheel  by 


p  =  1  +  7-  +  20  -  2 

K  2  r 


[7.18] 


where  9 (Figure  13a)  Is  related  to  slnkage  and  shearing  stress 
at  A  through  Equation  [6.34] 


9„  =  9  +  4>  =0  +  t!-  cos  1  t 

i  a  wc  a  2  w 


[7.19] 
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Summarizing,  the  stress  distribution  along  the  bow  portion 
of  the  wheel  is  given  by 

x  =  t  (AC,  -0  <  e  <  0)  [7.20] 

w  a 

a  =  -  p  -  sin  24>  -  -  1  +  ^  +  20  +  sin  24> 

w  wc  2  r  wc 

0  +  2<D 

+  2  — — g - —  0  (AC,-0g  <  0  <  0)  [7-21] 

a 

Equations  [7.20]  and  [7.21]  give  the  stress  distribution  in 

terms  of  0,  0  ,  0  and  x  . 

r  a  w 

The  last  point  to  be  clarified  is  that  of  the  rigid  wedge  at 
the  bow,  delimited  by  the  free-surface  and  the  slipline  AB 
(Figure  13a). 

The  pressure  at  A  is  (from  Equation  [7.17]  with  0  =  -0  ) 

3 

p  =  1  +  ^  +  20  +  20  +44>  [7.22] 

A  2  r  a  wc 

This  result  is  independent  of  the  assumption  of  linearity  of  the 
line  AC  in  the  characteristic  plane. 

Hill  (9)  has  shown  that  the  rigid  wedge  at  the  bow  is  ef¬ 
fectively  in  a  rigid  (elastic)  state  if 


FA  £  1  +  3-  - 


[7.23) 


HYDRONAUTICS ,  Incorporated 


-80- 

the  sign  equal  corresponding  to  a  plastic  state  given  by  a 
centered  fan  at  A. 

Hence.,  from  Equations  [7.22]  and  [7-23]  we  find  that 

0  <  -  2 0  -  30  [7.24] 

r  4  a  wc 

For  0^  larger  than  the  limit  of  Equation  [7.24]  the  soil  at 
bow  cannot  sustain  in  a  rigid  state  the  stresses  generated  by  the 
wheel.  Then  the  plastic  flow  extends  upstream  causing  probably 
bulldozing  effects. 

The  application  of  condition  [7.24]  to  the  relationship 

between  0  ,  0  and  <t>  in  Figure  20  shows  that  [7.24]  is  violated 
r  a  wc 

only  for  the  largest  sinkage  (0  =  45  ,  z  »  0.3)  and  at  small 

a 

shear  stress  between  the  wheel  and  the  soil  (low  slippage).  The 
corresponding  region  is  shaded  in  Figure  20.  The  conclusion  is 
that  the  picture  of  steady  flow  assumed  in  this  work  is  probably 
not  possible  for  driven  wheels  at  sinkages  larger  than  z  =  0.3. 

We  have  not  found  in  literature  [5L  [  1 8 ] ,  [22]  examples  of  tests 
in  which  larger  sinkage  values  have  been  realized. 
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8.  FORCES  ACTING  ON  THE  WHEEL  AND  THE 
MINIMUM  SLIPPAGE  COEFFICIENT 

Equations  [7.5],  [7.9],  [7.20]  and  [7. 21]  gave  the  approxi¬ 
mate  stress  distribution  on  the  wheel.  The  distribution  is 
represented  in  Figure  22  for  the  example  discussed  in  Section  6 
(0  =  30°,  z  =  0  134j  <J>  =  35°,  t  =0.34).  From  the  diagram 

of  Figure  20  it  is  found  that  the  rec  very  angle  is  0^  =  17.2  . 

The  normal  stress  is  discontinuous  at  C,  the  jump  due  to 

the  rigid  core  being  2  sin  2$wc.  As  the  shearing  stress  on  the 

wheel  increases  this  jump  diminishes.  The  normal  stress  c  has 

its  maximum  value  at  the  bow  and  drops  along  the  wheel.  The 

Isotropic  pressure  is  continuous  along  the  wheel.  For  the  same 

sinkage,  but  a  perfectly  rough  wheel  (r  =  l)  the  normal  stress 

w 

distribution  is  continuous  and  0 ^  is  larger,  0^  =  23  (Figure  22). 
The  shearing  stress  is,  by  assumption,  constant  along  AC  and  drops 
to  zero  from  C  to  R. 

The  forces  acting  on  the  wheel  may  be  easily  found  by 

integrating  the  stresses,  0  being  the  integration  variable  in 

the  range  -0  <  0  <  0  . 

a  r 

A .  The  Vertical  Force  W  (Flotation,  Figure  22) 

The  vertical  force  acting  on  the  wheel  has  four  components 


W  =  Wi  +  W2  +  W3  +  W4 


[8.1] 
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where  Wi  is  the  vertical  projection  of  the  normal  stress  a w 
acting  on  AC,  W2  results  from  on  CR,  W3  from  t  on  AC  and  W4 
from  t  on  CR.  All  these  components  but  W4  are  directed  upward 
(positive ) . 

Since  dx  =  cos  6  d0  and  dy  -  sin  0d0,  the  expressions  of  the 
W  components,  according  to  Equations  [ 7 • 5 1 >  [7 -9),  [7-20]  and 
[7.21]  are 


O  r 


Wi  =  -  I  = 


(1  +  f  +  2%  +  sin  2*wc.)  +  2 


0  0  +  <t> 
l  wc 

e.  -  <t> 

l  wc 


cos  Gde 


-0 


1  +  |  +  20  +  Sin  24>  j  sin  0  -  2  1  _™  —  (l-0asin0a-cos  0g ) 

^  r  w  I  £  wc 


[8.2] 


R  0 


W2  =  -  lawdx  =  (1  +  -  +  20r  -  sin  2$wc) 


\ 

i-f 

\  r  / 


cos  0d0 


1  -  cos  0 


=  (1  +  2  +  20r  '  Sin 


wc '  0 


[8.3] 
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•t  z  =  cos  2v  (1-cos  9  )  [8.4] 

w  wc  a 


Hence,  by  summation,  W  is  expressed  in  a  closed  form  as  a 

function  of  6  (sinkage),  $  or  t  (shearing  stress  on  the  wheel 

related  to  slippage),  6  (recovery  angle)  ar.d  9^.  Since  Q ^  is 

itself  a  function  of  9  and  t  (Figure  20)  and  the  same  is  true 

a  w  x 

for  9 .  (9 .  -  9  +  <t  )  the  flotation  is  in  fact  a  function  of  0 

l  v  l  a  wc'  a 

and  t  solely, 
w 

In  Figure  23  the  dependence  of  W  on  t  and  z  is  represented 

w 

graphically.  The  values  have  been  obtained  ty  using  Figure  20 
for  finding  9 and  Equations  [8  2] -[8b]  for  W. 

The  computations  show  that  the  major  part  of  VI  is,  of  course, 

given  by  VI  x  As  tne  shearing  stress  -v  and  the  sinkage  z  increase, 

w 

however,  the  other  components  of  W  'nay  reach  as  rn^cn  as  30  per¬ 
cent  of  its  value.  The  neglec^ion  of  the  normal  stresses  on  the 
rear  portion  as  assumed  by  Bekker,  theory  is  not  justified  in 
this  case. 
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As  expected  W  increases  sensibly  with  the  sinkage.  An  in¬ 
teresting  feature  is  the  drop  of  W,  for  a  given  sinkage,  with  t  , 
i.e.  with  the  slippage.  A  computation  of  an  "average"  stress  ob¬ 
tained  by  dividing  the  flotation  W  to  the  contact  length  6 +  0  , 
based  on  the  results  of  Figures  20  and  23,  shows  that  W/f3  +  0^) 
has  an  average  value  of  3*9  (for  15°  <  <  45°,  10°  <  <t>w  <  40  ) 

with  deviations  of  maximum  ±  13  percent. 

B.  Horizontal  Force  H  (Resistance,  Drawbar  Pull,  Figure  20) 

The  horizontal  force,  considered  positive  if  it  is  resistive 
and  negative,  if  it  is  propulsive,  has  again  four  components 


H  =  R  +  H2  +  H3  +  H4 


[8.6] 


where  R  results  from  a  on  AC,  H2  from  a  on  CR,  H3  from  t  on 

w  w 

AC  and  H4  from  T  on  CR  (Figure  22). 

The  expressions  of  the  different  components  of  H  are  as 
follows 


R  =  -jaw  dy  = 
A  < 


e.+  * 

u  W  C 

(1  +  tt/2  +  2 Br  +  sin  2<l>wc)  +  2  - - : -  6 


e 


<t> 


wc 


sin  0d0 


=  (1  +  tt/2  +  20p  +  sin  2$ wQ  )  ( 1 -cos  0Q  ) 


0  n  +  <t> 

- —  (0  cos  0  -  sin  0  ) 

0.  -  <t>  v  a  a  a 

l  vie 


2 


[8.7] 
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H2  =  - 


R 


a  dy  =  -  /  ( 1  +  tt/2  +  20  -  sin  2$  )  1  - 

w  J  I  '  '  r  wc '  0 


cos  0d0  = 


=  (l  +  tt/2  +  20  -  sin  20  ) 

v  '  r  wc' 


sin  0  \ 

r> 

-  -  1 


0 


[8.8] 


h3  = 


t  dx  =  -  t  /  cos  0d0  =  -  t  sin  0  =  -  cos  24>  sin  6 

w  w  |  w  a  vie  a 


-0 


[8.9] 


H4  =  - 


r 

/  xdx  =  - 

r  r 

0  \ 

T 

1  -  -T- 

J 

W 

0 

1  r  / 

cos  0d0  =  T 


cos  0  -  1 

_ r _ 

w  0 


=  cos  2<t> 


cos  0  -  1 

_ r _ 

wc  0 


[8.10] 


The  force  R  is  a  resistive  one,  the  other  three  forces  are 

propulsive.  Again,  for  a  given  sinkage  and  tw  the  values  of  H 

may  be  found  through  the  relationships  between  0^  and  0q,  tw. 

In  Figure  24  the  value  of  H  as  a  function  of  t  (or  <t>  )  for 

w  wc ' 

different  sinkages  is  represented,  as  well  as  the  value  of  R. 

In  usual  terminology  R  is  the  resistance,  H  is  the  minus  draw¬ 
bar-pull  and  H-R  is  the  traction. 
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As  expected  the  resistance  Increases  with  the  slnkage.  For 
a  given  slnkage,  both  total  forces  and  resistance  drop  as  the 
shear  stress  Increases,  but  the  total  force  drops  fasuer  becoming 
negative  (l.e.  drawbar  pull)  at  high  values  of  tw  and  low  values 
of  z . 


C.  The  Torque  M 

The  torque  acting  on  the  wheel  arises  only  from  shearing 
stress.  It  is  given,  therefore,  by 


D.  The  Minimum  Slippage 

The  relationship  between  the  assumed  constant  t  =  tw  along 
the  wheel  (on  AC)  and  the  slippage  has  been  discussed  In  Sec¬ 
tion  6d .  Only  a  detailed  analysis  of  the  boundary  layer  between 
the  wheel  and  the  plastic  zone  may  reveal  the  true  nature  of  this 
relationship.  At  this  moment  one  may  use  an  empirical  relation¬ 
ship  between  tw  and  S.  At  any  rate,  the  analysis  (Section  6d) 
has  shown  that  a  minimum  slippage  is  necessary  in  order  to  main¬ 
tain  positive  shear  stress  along  the  entire  wheel.  From  Equa¬ 
tion  [6.39]  S  is  given  by 
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Sin  e 


s 


m 


Sin 


wc 


4-  sin  6r 

C 


[8.12] 


i.e.  S  is  depending  on  sinkage  and  on  r  (or,  otherwise  stated 
m  w 

t  depends  on  and  sinkage). 


In  Figure  23  the  dependence  of  S  on  r  and  z  Is  repre- 

m  w 

sented  graphically.  For  a  given  sinkage  S  increases  with  t  , 

m  w 

the  increase  being  very  fast  at  high  t (?w  >  0.8) 

The  flow  pattern  suggested  in  this  report  for  a  driven  wheel 
is  possible  only  for  S  >  Sm  .  At  lower  slippage  coefficients 

the  wheel  starts  to  skid.  For  the  small  sinkages  encountered  in 

applications  and  moderate  z the  required  minimum  slippage  are  of 

the  order  0.35  -  0.40.  at  lower  values  the  wheel  i  no  more 

totally  slipping  and  a  different  solution  has  to  be  contemplated. 


For  driven  wheels  without  lateral  spuds  ,  a  slippage  coef¬ 
ficient  of  0.35  -  0.^0  is  not  high  and  it  may  be  said  that  the 
present  solution  covers  the  whole  practical  range  of  slippages 
for  driven  wheels  (l8). 
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9.  DISCUSSION  OF  RESULTS 

A .  Comparison  With  Experiments  (General) 

As  stated  in  Section  2b,  the  existing  experimental  results 
are  not  adequate  for  a  quantitative  verification  of  the  present 
work.  Adequate  experiments  for  this  purpose  have  to  fulfill  the 
following  requirements:  (i)  the  soil  behavior  has  to  be  con¬ 
trolled  and  known  (homogeneity,  stress-strain  relationship,  in¬ 
compressibility).  A  soil  behaving  as  close  as  possible  to  the 
plastic-rigid  model  has  to  be  used,  (ii)  the  wheel  shape  and 
motion  has  to  be  similar  to  those  assumed  here,  i.e.,  two-di¬ 
mensional,  steady  and  with  correct  scaling  of  the  inertial  num¬ 
ber  H  (ill)  accurate  measurements  of  velocity  distribution  in 
the  soil  body  and  stress  distribution  along  the  wheel  have  to  be 
carried  out. 

Since  no  experimental  data  satisfying  these  requirements 
have  been  available,  only  a  qualitative  comparison  between  the 
trends  shown  by  actual  experiments  and  the  theory  Is  possible. 
Quantitative  comparison  is  also  made  when  possible. 

B.  The  Recovery  Angle 

For  the  first  time  an  approximate  theoretical  method  of 
predicting  the  magnitude  of  the  recovery  angle  0 has  been  pre¬ 
sented.  The  values  of  0^  seems  high  when  compared  with  general 
accepted  values  of  0^  =  5°  -  6°.  One  has  to  keep  in  mind,  how¬ 
ever,  that  this  latter  value  is  based  on  experiments  with  wheels 
of  finite  width  in  which  case  a  rut  is  created  and  part  of  the 
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sol]  is  moved  sidewise.  Thus  causes,  obviously,  a  reduction  of 
0  .  Moreover,  if  the  soil  is  compressible  (as  assumed  in  Bekker's 
theory)  again  smaller  values  of  0 will  be  observed. 

In  the  experiment  of  (l8)  (Figure  25)  done  with  wide  rollers 
in  clay,  values  of  9 of  l6°  have  been  observed;  in  (22)  9 ^  seems 
to  be  even  higher,  so  that  our  results  are  by  no  means  beyond  the 
practical  range,  for  the  stipulated  condition. 

The  Increase  of  9 with  t^,  i.e.  with  slippage,  (Figure  20) 

seems  intuitively  correct.  One  would  expect  that  az  stall  0^ 

will  reach  its  maximum  value  0  =  0  .  In  our  results  0  is  below 

r  a  r 

0  even  for  S  —1.  It  must  be,  however,  emphasized  that  the 
a  m 

stall  cannot  be  regarded  as  a  limit  case  of  flow  and  experiments 
tend  to  confirm  it  (see  (21)). 

C.  The  Normal  Stress  Distribution  and  the  Flotation  W 

The  maximum  normal  stress  a  has  been  found  at  the  bow  A 

w 

(Figure  22  is  an  example).  The  measurements  of  (l8),  ( 1 5 )  show 
that  high  values  of  a exist  at  the  bow,  contradicting  Bekker- 
Bernsteln  assumption  (Figure  25).  The  maximum  is  close  to  the 
bow  in  (15),  hut  near  the  bottom  in  (l8)  Experiments  in  sand 
also  show  that  the  maximum  is  between  the  bottom  and  the  bow. 

In  discussing  the  work-hardening  effect  it  will  be  shown  that 

the  maximum  a  has  in  fact  to  be  displaced  from  the  bow,  but  it 
w 

will  be  closer  to  it  as  the  slippage  increases. 
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No  stress  discontinuity  has  been  observed  in  the  measure¬ 
ments  available  to  us,  but  a  gradual  drop  of  a from  the  bottom 
to  the  rear  has  been  measured  in  both  (15)  and  (l8)  (Figure  25). 

The  discontinuity  is  a  crude  approximation  for  this  high  drop. 

The  flotation  W  (Figure  23)  has  values  of  the  order  found  in 
(5b  being  somehow  higher.  Since  the  experiments  in  (5)  have 
been  carried  out  with  narrow  wheels,  the  discrepancy  is  expectable. 
An  effect  observed  in  the  same  experiments  (5)  but  ignored  by 
Bekker's  theory  and  empirically  considered  by  Reece  (13)  in  theory 
on  wheels  in  sand  is  the  drop  of  flotation  with  slippage.  The 
present  theory  has  predicted  this  effect  (Figure  23).  Moreover, 
its  reason  is  quite  transparent.  As  t^,  i.e.  the  slippage,  in¬ 
creases,  the  angle  8.-8  +  <t>  (Figure'  13a)  decreases  and  the 

)L  3  VJ  C 

slipline  AC  approaches  the  wheel.  Since  the  magnitude  of  the 

; 

pressure  at  the  bow  A  depends  on  the  bending  of  this  slipline, 
the  peak  value  of  a w  =  a decreases  as  the  slipline  approaches 
the  wheel.  The  above  effect  may  be  expressed  in  a  slightly  dif¬ 
ferent  form  by  asserting  that  for  a  given  W  the  sinkage  will  in¬ 
crease  with  slippage. 

The  normal  stress  vary  along  the  wheel,  which  is  different 

of  Uffelmann's  assumption  of  constant  cr  (l8).  The  average 

normal  stress  W/(0  +  8  ')  is  of  order  3.9  in  the  present  results 

3  r 

compared  with  the  value  of  order  5. 7  assumed  by  Uffelmann  (l8). 

Our  value  is  much  closer  to  his  measured  value  of  3.6  (Figure  25). 
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D.  Resistance  and  Drawbar  Pull 

The  high  normal  stress  at  rear  is  the  main  resistive  factor. 
Drawbar-pull  has  been  obtained  only  for  small  sinkages  and  high 
slippages  (Figure  24).  No  experimental  data  for  rims  have  been 
found,  excepting  (18). 

Experiments  with  towed  wheels  tend  to  show  that  irrespective 
of  the  type  of  tested  soil,  the  experimental  data  group  satis¬ 
factory  in  a  diagram  in  which  R/W  (resistance  over  flotation) 
are  plotted  as  a  function  of  sinkage.  Firth  (6)  has  found  from 
different  measurements  for  cylindrical  wheels  that  a  good  fit 
for  towed  rigid  wheels  is 


[9.1] 


with  K  =  0.75  and  n  =  0.45. 

Based  on  the  results  of  Figures  23  and  25  the  values  of  R/W 

as  a  function  of  z  for  different  t  have  been  plotted  on  a 

w 

logarithmic  paper  (Figure  26).  The  dependence  is  well  approximated 

by  straight  line  which  have  approximately  the  same  slope  of  n=0.5. 

The  value  of  K  differs  with  slippage,  being  lower  at  high  slippages. 

For  t  =  0.19*  i.e.  relatively  low  shear  stress,  K  =  0.71.  Hence, 
w 

the  present  theory  provides  values  quite  close  to  those  found  for 


towed  wheels. 
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E.  Influence  of  Work-Hardening,  Compressibility  and  Side  Effects 

It  is  a  well  established  fact  that  clay  shows  a  dependence 
of  the  yield  stress  on  the  strain  which  resembles  work-hardening, 
i.e.  an  increase  of  the  yield-stress  with  the  strain. 

The  qualitative  influence  of  the  work-hardening  has  been  con¬ 
sidered  in  different  sections.  It  is  worthwhile  to  recall  these 
effects  and  to  discuss  them  briefly  here. 

A  first  effect  is  the  diffusion  of  the  discontinuity  line  AC 
which  becomes  a  zone  of  transition  rather  than  a  line.  At  this 
stage  it  is  difficult  to  estimate  the  thickness  of  the  transition 
layer,  but  this  may  explain  some  upstream  deformations  not  en¬ 
countered  in  a  rigid-plastic  model. 

Another  effect  is  reflected  in  the  normal  stress  distribu¬ 
tion.  Since  the  pressure  increases  with  strain,  it  is  expected 
that  at  the  bow  A  the  pressure  will  be  lower  than  the  value  given 
by  the  rigid-plastic  model  and  will  increase  rapidly  along  the 
wheel,  reaching  the  rigid-plastic  value  in  the  region  of  high 
strain. 

The  compressibility  of  the  soil  will  cause  changes  in  the 
flow  pattern.  An  apparent  effect  is  the  creation  of  a  rut  (i.e. 
the  diminution  of  0  )  and  additional  dissipation  in  compaction. 

A  narrow  wheel  creates  a  three-dimensional  flow  pattern. 
Experiments  (21)  show  that  the  flow  beneath  the  wheel  exhibits 
the  same  features  as  in  the  case  of  a  roller.  A  rut  is  now, 
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however,  possible.  For  this  reason  and  due  to  side  effects  the 
flotation  will  probably  be  lower  than  in  the  two-dimensional  case. 

Future  theoretical  work  which  will  incorporate  these  effects 
is  suggested  in  section  11. 

F.  Energy  Considerations 

The  equation  of  energy  is  in  the  case  of  an  incompressible 
material  a  simple  transcription  of  the  momentum  equations 
and  it  is,  therefore,  not  very  helpful  in  the  study  of  the  me¬ 
chanics  of  motion.  It  is  useful  for  finding,  for  instance,  the 
temperature  changes  in  the  material-which  involves  additional 
empirical  data  on  thermal  properties. 

Since  it  is  usual  to  express  some  results  found  from  the 
momentum  equations  in  terms  of  energy,  it  will  be  worth¬ 
while  to  examine  briefly  the  energy  balance  in  the  case  of  soil- 
wheel  interaction.  Considering  the  soil  body  (Figure  l),  one  can 
apply  the  equation  of  mechanical  energy  to  a  control  volume 
delimited  by:  a  vertical  section  1-1  far  upstream,  a  vertical 
section  2-2  far  downstream  (  both  from  the  free-surface  to 
infinity  downward),  the  free-surface  and  the  soil-wheel  inter¬ 
face.  One  of  the  basic  assumptions  underlying  this  work  is  that 
of  uniform  flow  and  unperturbed  free-surface  far  from  the  wheel. 

In  this  case  the  energy  flux  through  sections  1-1  and  2-2  is 
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d“lt  "  V  !  (oxl’  -  °x2,)dy'  [9'Sl 

The  free-surface  is  stress-free;  there  is  no  energy  flow 
across  it. 

Work  is  done  on  the  soil-wheel  interface  by  shearing  stress 
solely.  The  work  rate  is 


dE2 

dt 


V  '  /  Tds 
w  J 


The  energy  equation  states  that 


dEi  dE2  _  dEd 

dt  dt  "  dt 


[9.3] 


[9.4] 


where 


dEd 

dt 


is  the  rate  of  energy  dissipation. 


From  the  momentum  equation 


f(cx i'  -  ax2-  )dy *  =  H'  [9.5] 

where  H'  is  the  total  horizontal  force  acting  on  the  wheel.  By 
definition 


/  Tds  =  -  M 1 


[9.6] 
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where  M'  is  the  torque  acting  on  the  wheel.  Also  V'  =  U)'r'. 

V) 

Hence,  Equation  [9.^]  becomes 


H'U  1  -  M'co'r'  =  [9.7] 

o  dt 

Since  in  the  case  of  a  driven  wheel  the  torque  M'  is  applied 
to  the  wheel  in  clockwise  direction,  M'  <  0.  H'  is  positive  if  it 
is  a  resistance  and  negative  if  it  is  a  drawbar  pull. 

Equation  [9. 7]  states  the  simple  result  that  the  work  done 
by  the  horizontal  force  and  the  torque  are  equal  to  the  dissipa¬ 
tion  . 

Dissipation  occurs  in  the  plastic-rigid  model  in  the  plastic 
zone  solely.  This  zone  has  to  be  divided  in  three  regions: 
the  narrow  transition  layer  between  the  rigid  and  plastic  zone 
(considered  a  discontinuity  layer},  the  bulk  of  the  plastic  zone 
and  the  boundary  layer  at  the  soil  wheel  interface. 

The  dissipation  in  the  transition  layer  //  t  '  (bV/dn  )dnds 
where  s  and  n  are  directions  along  and  normal  to  the  transition 
layer  respectively  and  V'  is  the  tangential  velocity.  For  a 
perfect  plastic  material  and  no  inertial  effects  t'  =  k'  and 
f  (bV1  /bn)  dn  =  6V'  (the  velocity  jump).  The  matter  is  more  com¬ 
plex  when  work  hardening  and  inertial  effects  are  considered, 
but  it  will  not  be  discussed  here. 
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The  dissipation  in  the  plastic  zone  has  been  discussed  in 
Section  3F* 

In  the  boundary  layer  near  the  wheel  the  dissipation  equals 
/ T'(o)'r'  -  V')  ds,  where  V'  is  the  soil  velocity. 

In  any  computations  based  on  energy  considerations  special 
attention  has  to  be  paid  to  the  dissipation  in  the  two  thin  layers 
which  may  account  for  a  large  part  of  the  total  dissipation. 

G.  Extension  of  the  Method  to  Different  Stress  Conditions 
Along  the  Wheel 

The  solution  of  the  plastic  flow  in  the  bow  region  (Sec¬ 
tion  6d)  has  been  obtained  by  using  the  condition  of  constant 
shear  stress  t  =  tw  along  the  wheel  on  AC. 

There  is  no  major  difficulty  in  solving  the  problem  with 
an  arbitrary  distribution  of  t  along  the  wheel,  provided  that 
the  same  picture  of  the  slipline  field  is  adopted. 

In  section  6d  the  t  boundary  condition  was  used  in  order 
to  locate  the  image  of  the  line  AC  in  the  characteristic  plane 
(Figure  14). 

Assume  that  a  relationship  between  t  and  0  (position  on  the 

wheel  -6  <  0  <  0)  is  given 

3 


t  =  r(e) 


[9.8] 
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Equation  [4.7]  relates  t  to  t  and  0.  Hence,  from  Equa¬ 
tions  [4.7]  and  [9*8]  we  have 

x(0)  =  cos  2  ($-0)  [9*9] 


which  permits  to  determine  $-0  along  the  wheel  as  a  function 
of  <t> .  From  the  other  hand  Equations  [4.8]  which  reflects  the 
fact  that  the  wheel  is  a  streamline,  gives 


v 

=  -  tg  (<t>-P) 


a 


[9.10] 


Since  <t>  is  related  to  a,  /3  by  Equation  [6.40],  v  /v  on  the 

O  OL 

wheel  may  be  computed  for  a  given  x(0).  The  same  procedure  as 
in  the  case  of  constant  t  permits  the  location  of  the  AC  line  in 
the  characteristic  plane  (Figure  14 ) . 


This  general  procedure  applies  to  frictional  conditions,  for 
instance,  in  which  x  is  proportional  to  a  .  Since  is  a  known 
function  of  a,  /3  agaln?  is  relatively  easy  to  find  the  points  in 
the  a,j3  plane  which  satisfy  the  assumed  condition. 


A  dependence  of  x  along  the  wheel  on  the  strain  of  the 

thin  plastic  layer  adjacent  to  the  wheel  may  also  be  treated  in 

the  same  way.  The  strain  is  depending  on  the  difference  between 

the  wheel  velocity  V  and  the  plastic  velocity,  and  the  location  0. 

w 

If  x  is  a  known  (say  empirical)  function  of  strain,  again  a  re¬ 
lationship  between  v^  and  may  be  formed  and  the  position  of 
AC  may  be  determined. 
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Concluding,  as  long  as  the  stress  condition  along  the  wheel 

may  be  translated  in  a  kinematic  condition  for  v  ,  v„,  ♦  -  it  is 

a  /3 

relatively  easy  to  solve  the  plastic  flow  problem  in  the  bow 
region . 

At  the  present  stage,  because  of  the  lack  of  sufficient  ex¬ 
perimental  or  theoretical  evidence,  it  seems  unjustified  to  go 
beyond  the  assumption  of  constant  r . 
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10.  SUMMARY  AND  CONCLUSIONS 

A  comprehensive  theoretical  solution  of  the  soil-wheel 
interaction  is  for  the  first  time  presented.  The  soil  flow  is 
steady  and  two-dimensional,  the  wheel  is  cylindrical  and  rigid. 
Only  the  case  of  a  driven  slipping  wheel  is  considered.  The 
soil  is  assumed  to  be  incompressible  and  to  behavellkea  rigid- 
plastic  material,  which  may  be  approximately  valid  in  the  case 
of  a  soft  saturated  clay. 

The  flow  of  the  soil  in  the  region  beneath  the  wheel  is 
solved  by  integrating  the  equations  of  motion  with  appropriate 
boundary  conditions.  For  the  first  time  the  influence  of  the 
Inertial  terms  is  discussed.  Although  the  solution  is  obtained 
with  the  aid  of  the  quasi-static  equations,  the  flow  field  is 
built  such  that  inertial  influence  is  taken  into  account  at 
discontinuity  lines  and  points  where  the  Inertial  effects  are 
the  most  important. 

In  essence  a  plastic  soil  developes  high  normal  stresses 
on  the  wheel,  responsible  for  flotation  and  resistance,  in  the 
following  way:  The  soil  flowing  beneath  the  wheel  is  under¬ 
taking  large  deformations  which  cause  a  plastic  flow;  in  a  plasti 
region  the  maximum  shear  stress  is  everywhere  equal  to  the  yield 
stress;  the  equilibrium  of  an  element  moving  along  curved  stream¬ 
lines  requires  pressure  gradients  in  order  to  sustain  the  shear 
stress.  These  pressure  gradients  build  up  in  high  pressures  be¬ 


neath  the  whee 
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An  approximate  solution  of  the  plastic  flow  problem  is  pre¬ 
sented.  The  recovery  angle,  the  vertical  and  horizontal  forces 
acting  on  the  wheel  are  found  in  an  analytical  form  as  functions 
of  the  sinkage  and  the  magnitude  of  the  shear  stress  on  the 
wheel  boundary.  The  latter  is  related  to  the  minimum  slippage 
coefficient,  necessary  to  maintain  positive  shear  stress  on  the 
wheel . 

The  solution  permits  for  given  wheel  radius,  vertical 
load  and  shear  stress  on  the  wheel  (or  minimum  slippage)  to  de¬ 
termine  the  sinkage  and  the  horizontal  force.  The  procedures, 
using  the  analytical  expressions  or  the  graphs,  is  the  following 
for  given  vertical  load  (W1  )  radius  (r1)  yield  stress  (k1 )  and 
shear  stress  along  the  wheel  (t'w)  (or  minimum  slippage  )  the 
sinkage  z=z'/r'  found  from  Figure  23.  Figure  24  permits  the 
determination  of  the  resistance  R'  and  the  total  horizontal 
force  HI  .Figure  20  gives  the  recovery  angle  9^. 

The  results  exhibit  all  the  trends  well  known  from  experi¬ 
ments  (like  dependence  of  flotation  on  sinkage  and  slippage), 
without  using  the  artificial  empirical  assumptions  adopted  in 
previous  theories. 

Although  some  quantitative  agreement  between  the  present 
theoretical  results  and  measurements  has  been  found,  more  and 
better  experiments  are  necessary  in  order  to  validate  the  theory 
For  this  reason  the  present  work  has  to  be  regarded  as  a  pre¬ 
liminary  study. 
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The  plastic-rigid  theory  has  proved  itself  as  an  useful 
tool  of  analysis  of  soil-wheel  Interaction.  The  results  so  far 
obtained  encourage  the  continuation  of  theoretical  investiga¬ 
tions  which  will  further  advance  the  understanding  of  the  vehicle 
dynamics  and  provide  sound  foundations  to  a  field  of  applied  me¬ 
chanics  dominated  by  empiricism  in  the  past. 
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11.  SUGGESTIONS  FOR  FUTURE  INVESTIGATIONS 

The  review  of  the  literature  and  the  theoretical  analysis 
carried  out  in  this  report  suggest  some  directions  of  special 
interest  to  be  followed  in  future  investigations  of  soil-wheel 
interaction  in  soft  soil. 

In  a  broad  sense  it  seems  that  three  directions  are  of  major 
importance:  (i)  work  toward  a  better  quantitative  characteriza¬ 

tion  of  different  soils  important  in  mobility  research  A 
mechanical  characterization  means  mainly  constitutive  equations 
relating  the  stress  and  strain  tensors  including  strain-rate  ef¬ 
fects.  The  rigid-plastic  equations  represent  one  of  the  simplest 
types  of  models,  (ii)  thorough  laboratory  work  with  wheels  moving 
in  carefully  prepared  and  controlled  soil  with  accurate  measure¬ 
ments  of  velocity  field  and  stresses  on  the  wheel.  Experimental 
study  of  the  boundary  conditions  at  the  soil-wheel  interface. 

(iii)  theoretical  work  in  which  the  soil-wheel  problem  is  attacked 
from  basic  principles,  using  existing  knowledge.  This  type  of 
work  has  to  incorporate  new  results  achieved  in  directions  (i) 
and  (ii). 

The  work  presented  here  pertains  to  point  (iii),  In  a  narrow 
sense,  the  following  theoretical  aspects  of  soil-wheel  interaction 
direct  continuations  of  the  present  work  may  be  contemplated 

^1)  Refinement  of  the  solution  in  order  to  solve 
analytically  the  flow  in  the  whole  plastic  region. 
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(2)  Study  of  the  Inertial  effects  by  solving  the  second 
order  approximation  In  the  perturbation  scheme. 

(3)  Investigation  of  work-hardening  effects. 

(4)  Extension  of  solution  to  towed  wheels. 

(5)  Study  of  motion  In  sand  soil. 

(6)  Study  of  three-dimensional  effects  for  wheels 
of  finite  width. 

(7)  Study  of  non-steady  effects  by  incorporating 
local  accelerations  and  the  time  factor  in  the  equations  of  flow. 

Because  of  the  complexity  of  the  wheel  problem  we  feel  that 
the  best  way  to  estimate  the  above  effects  will  be  a  study  of  a 
tracked  vehicle,  which  Is  geometrically  simpler.  We  strongly 
recommend  therefore,  that  future  work  should  stari  with  the  so¬ 
lution  of  plastic-rigid  flow  beneath  a  track  with  consideration 
of  Inertial,  work-hardening,  three-dimension^ i  and  unsteady  ef¬ 
fects.  The  next  immediate  step  should  be  the  extension  to  the 
wheel  problem. 
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APPENDIX  1 


THE  VELOCITY  COMPONENTS  AND 
0IN  THE  FREE -SURFACE  REGION 
(FIGURE  12)  AS  FUNCTION  OF  0,11 


(RAD.)  28--  v.  V  ♦(DEGREES)  .(RAD.)  2B--  *«  V  ♦(DEGREES) 
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APPENDIX  2 


THE  VELOCITY  COMPONENTS  AND  <P  IN 
THE  BOW  PLASTIC  REGION  (FIGURE  12) 
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APPENDIX  3 


THE  MAPPING  OF  THE  BOW  PLASTIC  ZONE 
(FIGURE  12)  FROM  THE  CHARACTERISTIC  PLANE 
ON  THE  PHYSICAL  PLANE 
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a.  IDEAL  ELASTIC-PLASTIC  MATERIAL  b. 
(FOR  RIGID-PLASTIC  y  '  =  0  ) 


T' 


WORK-HARDENING  MATERIAL. 

(FOR  RIGID  WORK -HARDENING  Y '  =0, 

FOR  RIGID-PLASTIC  y'e  =  0,  y'h  =  0) 


FIGURE  2  -  TH1  DEPENDENCE  OF  THE  MAXIMUM  SHEAR  STRESS  ON  STRAIN 
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FIGURE  3  -  STRESS  AND  VELOCITY  COMPONENTS  IN  CARTESIAN  COORDINATES 


a.  CARTESIAN 
COORDINATES 


b.  MAXIMUM-SHEAR 
STRESS  DIRECTIONS 


c.  PRINCIPAL 
DIRECTIONS 


FIGURE  4  -  STRESSES  AT  A  POINT  REFERRED  TO  DIFFERENT  DIRECTIONS 
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a.  COMPRESSION  b.  TENSION 

FIGURE  5  -  PLASTIC  FLOW  ALONG  A  FREE  SURFACE 


x 

a.  SLIPPING  BOUNDARY 
(V  >Vw) ,  RETARDED  FLOW 

(!r  <  0)'  n/A>*  *°>0 


b.  SLIPPING  BOUNDARY  (V  >V  ) 

w 

FLOW  >0),  o  >*-e>  - 


c.  SKIDDING  BOUNDARY  (V  <V  ) 

w 


0  -  0 
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FIGURE  6  -  PLASTIC  FLOW  ALONG  A  RIGID  BODY 
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tangential  velocity 
distribution 


a.  THIN  PLASTIC  LAVER 


a 


6V'=  v'  -  U  '  cos  <p 
a  o 

b.  KINEMATICS  OF 
DISCONTINUITY 


c.  SHEAR  STRESS 


FIGURE  7  - 


A 

A 


UNE  OF  VELOCITY  DISCONTINUITY  AT  THE  BOUNDARY  BETWEEN 

rigid  and  a  plastic  region 
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b.  A  CENTERED  FAN  AT 
A  CONCAVE  CORNER 


c.  FLOW  DETAIL  IN  SHEET-DRAWING  d.  FLOW  AS  IN  c.  ,  WITH  INERTIAL 

(FROM  [8  p.  173]).  AD-  DISCON-  TERMS  CONSIDERED.  (AD-DIS- 

TINUITY  LINE,  DAC  -  CENTERED  FAN  CONTINUITY  LINE ) 


FIGURE  8  -  FLOW  IN  A  CENTERED  FAN  AND  AT  THE  INTERSECTION 
BETWEEN  THE  FREE-SURFACE  AND  A  RIGID  BODY 
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plastk  zone  1 


plastic  zone  2 


FIGURE  9  -  A  LINE  OF  STRESS  DISCONTINUITY 
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'  shear  stress  acting  on 
^  the  rigid  material 


b.  SUP-LINE  FIELD 


FIGURE  10  -  QUALITATIVE  REPRESENTATION  OF 
THE  PLASTIC  FLOW  PATTERN 
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the  rigid  material 


FIGURE  11  -  AN  INCORRECT  SLIP-LINE  FIELD 


b.  THE  CHARACTERISTIC  PLANE  a,li 

FIGURE  12  -  THE  FREE-SURFACE  REGION 
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a.  THE  PHYSICAL  PLANE  x,y 


b.  THE  VELOCITY  TRIANGLE  AT  A 


c.  THE  CHARACTERISTIC  PLANE  a,H 


FIGURE  13  -  PLASTIC  FLOW  IN  THE  BOW  REGION 


-0.8  -0.6  -0.4 
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FIGURE  15  -  THE  VELOCITY  DISTRIBUTION  ALONG  AC  (  SOIL-WHEEL  INTERFACE) 
AND  AB  (MARGINAL  CHARACTERISTICS) 


8’0- 


FIGURE  16  -  THE  SUP-LINE  FIELD  IN  THE  BOW  REGION  (EXAMPLE  OF  EQUATION  6.53) 
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b.  CHARACTERISTIC  PUNE 


FIGURE  19  -  CONTINUITY  OF  FLOW  IN  REGION  CDERC 
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vW  -DEGREES 


FIGURE  20  -  THE  DETACHMENT  ANGLE  B  AS  FUNCTION  OF  SINKAGE 

r 

AND  SHEAR  STRESS  ON  THE  WHEEL 
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FIGURE  21  -  THE  EQUILIBRIUM  OF  THE  RIGID  CORE  CBDC 


FIGURE  22  -  THE  APPROXIMATE  DISTRIBUTION  OF  STRESS  ALONG 
THE  WHEEL  (EXAMPLE  OF  EQUATION  6.53) 


FIGURE  24  -  RESISTANCE  AND  TOTAL  HORIZONTAL  FORCE  AS  FUNCTIONS 
OF  SHEAR  STRESS  AND  SINKAGE 
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DIRECTION  OF  TRAVEL 


W’  6849  x  180 


k'  r'  (9  +  6  )  6.61  x  27  x  12  x  61  x  3.14 

a  r 

(b) 

FIGURE  25  -  MEASURED  NORMAL  STRESS  DISTRIBUTION 
(a)  REF.  (15), p.  27  (b)REF.  (18),  p.  120 
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FIGURE  26  -  RESISTANCE  OVER  FLOTATION  AS  FUNCTION  OF  SINKAGE 
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